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Abstract 

Concentrated  photovoltaic  (PV)  technology  represents  a  growing  market  in  the 
field  of  terrestrial  solar  energy  production.  As  the  demand  for  renewable  energy 
technologies  increases,  further  importance  is  placed  upon  the  modeling,  design,  and 
simulation  of  these  systems.  Given  the  U.S.  Air  Force  cultural  shift  towards  energy 
awareness  and  conservation,  several  concentrated  PV  systems  have  been  installed  on  Air 
Force  installations  across  the  country.  However,  there  has  been  a  dearth  of  research 
within  the  Air  Force  devoted  to  understanding  these  systems  in  order  to  possibly  improve 
the  existing  technologies.  This  research  presents  a  new  model  for  a  simple  concentrated 
PV  system.  This  model  accurately  detennines  the  steady  state  operating  temperature  as  a 
function  of  the  concentration  factor  for  the  optical  part  of  the  concentrated  PV  system,  in 
order  to  calculate  the  optimum  concentration  that  maximizes  power  output  and 
efficiency.  The  dynamic  thennal  model  derived  is  validated  experimentally  using  a 
commercial  polysilicon  solar  cell,  and  is  shown  to  accurately  predict  the  steady  state 
temperature  and  facilitates  computer  analysis  and  prediction  of  the  ideal  concentration 
factor. 
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THERMAL  MODEL  OF  A  CONCENTRATED  PHOTOVOLTAIC  SYSTEM 


I:  Introduction 

Energy  is  the  key  enabler  to  successful  military  operations  on  a  global  scale.  To 
fly,  fight  and  win  in  air,  space,  and  cyberspace  requires  the  U.S.  Air  Force  to  use  vast 
amounts  of  energy.  Energy  is  the  fundamental  requirement  to  execute  every  Air  Force 
mission  and  the  role  that  energy  plays  in  successful  military  operations  from  a  tactical, 
strategic,  and  operational  perspective  cannot  be  understated.  As  the  largest  consumer  of 
energy  in  the  Department  of  Defense,  energy  awareness  within  the  Air  Force  is 
paramount  to  continued  successful  operations.  Consequently,  energy  issues  regarding  the 
security  and  sustainable  availability  of  this  resource  have  created  an  especially  acute  need 
to  evaluate  energy  resources  in  the  context  of  military  operations. 

Increasingly,  renewable  sources  of  energy  such  as  wind,  solar,  and  hydroelectric 
are  being  evaluated  for  their  applicability  to  certain  military  energy  system  requirements. 
Massive  solar  power  systems  have  been  installed  on  Air  Force  installations  across  the 
country  to  supplement  on-base  power  demands.  Furthermore,  alternative  fuels  are  being 
developed  and  tested  in  Air  Force  laboratories  to  diminish  dependence  upon  non¬ 
renewable  energy  sources,  so  it  is  evident  that  energy  issues  within  the  Air  Force  prove 
vitally  important  [1].  However,  in  order  to  obtain  a  well-rounded  objective  perspective 
on  energy  issues  within  the  Air  Force,  it  is  useful  to  evaluate  global  energy  issues. 
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1.1  Global  Energy  Issues 

Energy  is  a  primary  catalyst  to  development  and  globalization  in  the  modem 
world.  Furthermore,  as  the  global  population  continues  to  increase  the  demand  for 
energy  will  continue  to  rise  an  estimated  44  percent  over  the  current  level  by  2030  [1].  A 
growing  sector  of  energy  production  is  based  upon  renewable  energy  sources  such  as 
wind,  hydroelectric,  and  solar  power.  The  fundamental  problem  with  energy  production 
today  is  that  the  majority  of  energy  sources  have  a  finite  long-tenn  supply.  In  the  United 
States  during  2009,  according  to  the  U.S.  Energy  Information  Administration,  nearly  all 
of  the  electricity  consumed  was  produced  from  the  non-renewable  fossil  fuels  oil,  coal, 
and  natural  gas,  shown  below  in  Figure  1. 


U.S.  Energy  Consumption  by  Energy  Source,  2009 


Total  =  94.578  Quadrillion  Btu 


Total  =  7.744  Quadrillion  Btu 


--Solar  1% 

—  Geothermal  5% 
-Biomass  waste  6% 

Wind  9% 

Biofuels  20% 


Wood  24% 


Hydropower  35% 


Figure  1:  Global  Energy  Consumption  [2] 

As  the  global  population  continues  to  deplete  fossil  fuels  from  the  earth,  an 
imminent  need  is  placed  upon  the  exploitation  and  production  of  alternative  renewable 
sources  of  energy. 

The  solar  energy  resource  is  one  of  the  most  promising  solutions  to  the  global 
energy  supply  problem.  Every  day  the  sun  radiates  to  the  earth  roughly  10,000  times  the 
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energy  used  by  the  global  population  [3].  This  massive  amount  of  energy  could 
potentially  make  fossil  fuel  consumption  obsolete  in  the  future,  but  further  research  and 
development  of  solar  energy  systems  is  necessary  to  make  the  cost  of  solar  energy 
competitive  with  conventional  non-renewable  sources  of  energy. 

1.2  Air  Force  Energy  Issues 

The  2010  Air  Force  Energy  Plan  has  set  forth  a  series  of  energy  related  goals, 
with  the  primary  theme  “ make  energy  a  consideration  in  all  that  we  do ”  [4],  Figure  2 
shows  the  Air  Force  energy  consumption  relative  to  the  other  military  service  branches. 

U.S.  Armed  Forces  Fuel  Utilization 

(percent  of  total  fuel  cost) 

■  Air  Force 
|  Army 
|  Navy/Marines 

Figure  2:  Military  Energy  Consumption  [4] 

The  Air  Force  Energy  Plan  recognizes  energy  as  an  operational  enabler  and  also 
demonstrates  the  opportunity  for  Air  Force  centric  research  and  development  of  these 
systems.  Air  Force  infrastructure  planning,  specifically,  has  cited  the  need  for  renewable 
energy  systems  to  supplement  installation  energy  requirements  across  the  country  [5]. 
Furthermore,  at  deployed  and  forward  operating  locations,  where  the  raw  cost  of  fuel 
does  not  convey  the  cost  in  manpower  and  logistics  necessary  to  transport  this  resource  in 
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hostile  areas,  the  potential  implementation  of  modular  self-sustaining  solar  energy 
systems  to  power  communications  and  surveillance  equipment  proves  attractive. 
Engineers  working  on  the  Military  Electric  Project  (MEP),  a  joint  project  between  the 
U.S.  Army  and  National  Renewable  Energy  Lab,  are  currently  developing  a  rapidly 
deployable  solar  energy  system  specifically  for  troops  in  forward  operating  locations 
where  conventional  fuel  use  is  not  appropriate  [6]. 

The  idea  of  energy  sufficiency  for  troops  in  combat  reaches  far  deeper  than  the 
technological  details  indicate.  Commanders  in  Iraq  and  Afghanistan  have  also 
recognized  the  potential  lifesaving  impact  that  implementing  renewable  energy  systems 
into  combat  operations  could  have.  The  conventional  fuel  used  for  generating  electric 
power  in  combat  is  diesel  fuel,  which  drives  the  generators  powering  most  of  the 
deployed  military  facilities.  Since  the  most  dangerous  missions  today  in  tenns  of 
mortality  are  fuel  convoy  missions,  decreasing  the  number  of  fuel  convoys  necessary  to 
meet  the  enonnous  energy  demand  would  most  certainly  contribute  to  saving  lives  on  the 
battlefield  [7].  This  research  is  motivated  by  the  explicit  need  for  solar  energy  systems  in 
the  context  of  Air  Force  operations,  and  covers  the  theoretical  concepts  of  solar  cell 
operation  building  towards  a  new  solar  cell  model. 

1.3  Solar  Power  Basics 

The  theory  behind  solar  cell  operation  will  be  discussed  comprehensively  in  the 
body  of  this  research,  but  it  is  important  to  have  an  initial  top-level  understanding  of  solar 
energy  systems.  Solar  cells  generate  electricity  from  the  sun.  Photons  of  light  are 
absorbed  to  photogenerate  charge  carriers  within  the  cell.  The  term  “photovoltaic”  refers 
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to  the  photovoltaic  effect  discovered  by  Alexandre-Edomond  Becquerel,  an  experiment 
in  which  photons  of  light  falling  upon  a  material  generate  a  voltage  and  produce  a 
resulting  current,  and  often  times  solar  power  systems  are  referred  to  as  photovoltaic 
(PV)  systems.  While  not  every  solar  energy  system  uses  the  photovoltaic  effect  for 
energy  production  (such  as  solar  thennal  systems),  this  analysis  will  solely  cover 
photovoltaic  energy  systems. 

To  generate  usable  amounts  of  power,  individual  solar  cells  are  often  connected  in 
series  and  parallel  to  form  solar  arrays  that  produce  power  to  an  electrical  load. 
Applications  for  PV  systems  range  from  small  consumer  electronics  to  large-scale 
industrial  power  plants  generating  several  megawatts  of  power.  Solar  cells  are  also  used 
in  space-based  applications  to  generate  solar  power  on  satellites  and  other  space  vehicles. 
The  practical  modern  applications  for  PV  systems  are  obvious,  but  as  research  continues 
and  the  technology  is  developed  further,  solar  energy  as  a  primary,  sustainable,  and 
reliable  large-scale  energy  supply  may  become  reality. 

1.4  Concentrated  Solar  Power 

A  growing  field  within  the  photovoltaic  industry  is  concentrated  solar  power 
(CSP).  In  CSP  systems,  parabolic  mirrors  or  non-imaging  plastic  lenses  focus  sunlight 
onto  the  solar  cell.  The  principle  of  CSP  systems  is  that  the  focusing  optical  components 
increase  the  intensity  of  the  sunlight  on  the  cell  relative  to  the  degree  of  the  focus  factor 
of  the  optics,  increasing  the  photocurrent  through  the  cell  and  hence  the  overall  power 
generated  by  the  system.  The  immediate  advantage  of  CSP  systems  is  that  relatively 
inexpensive  optical  components  (mirrors,  lenses)  can  be  used  in  place  of  equivalent 
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amounts  of  expensive  semiconductor  material,  thus  theoretically  reducing  the  overall  cost 
of  these  systems. 

Most  modern  CSP  systems  use  an  optical  concentration  ratio  of  up  to  multiple 
hundreds  of  suns  [8].  Since  the  amount  of  semiconductor  material  necessary  is  relatively 
small  compared  to  the  physical  surface  area  of  the  optics,  high-efficiency  (but  expensive) 
semiconductor  material  can  be  used  for  photovoltaic  conversion,  and  CSP  systems 
represent  the  highest-efficiency  solar  energy  systems  currently  on  the  market  [8], 

The  increased  solar  concentration  creates  an  enormous  amount  of  heat  within  the  cell  so 
an  active  cooling  system  must  be  incorporated  into  every  high-concentration  ratio  system. 
Additionally,  since  the  focusing  optics  require  that  the  sun  radiate  at  a  90°  incidence 
angle  for  maximum  absorption,  a  solar  tracking  system  is  required  to  track  the  sun  across 
the  horizon  during  the  day.  The  tracking  system  further  increases  the  mechanical 
complexity  of  these  systems. 

Despite  the  mechanical  complexity  and  significant  heat  generated  in  CSP 
systems,  semiconductor  growth  and  processing  improvements  have  catalyzed  the  wide 
scale  implementation  of  CSP  arrays  for  terrestrial  energy  production.  Therefore,  it  is  of 
great  benefit  to  the  Air  Force  and  scientific  community  at  large  to  better  understand  these 
systems  at  the  theoretical  level  in  order  to  progress  this  technology  towards  meeting  the 
goal  of  military  energy  security  and  independence  and  decreasing  dependence  upon 
conventional  non-renewable  energy  sources. 

1.5  Research  Focus 

The  purpose  of  this  research  is  to  improve  the  means  for  modeling  and  to  develop 
a  method  for  designing  concentrated  photovoltaic  systems.  The  current  dearth  of 
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analytical  research  behind  these  optical  concentrator  systems  has  presented  an  exciting 
opportunity.  Instead  of  relying  upon  empirical  data  to  determine  the  optimum 
concentration  factor,  the  overall  goal  of  this  research  is  to  develop  a  design  tool  that  can 
accurately  calculate  the  ideal  optical  concentration  factor  based  on  extensive  dynamic 
thennal  modeling  of  the  device. 

The  foundation  of  this  research  lies  in  solid  state  semiconductor  physics. 

However,  extensive  research  has  also  been  conducted  in  the  areas  of  photonics,  thennal 
dynamics,  and  mathematics.  Furthermore,  MATLAB  is  used  to  generate  the  complex 
multi-order  system  models  to  accomplish  the  research  goal.  In  particular,  the  following 
issues  and  questions  are  addressed: 

1 .  What  is  the  mechanism  that  a  solar  energy  system  uses  to  convert  light  into 
electricity? 

2.  How  do  the  physical  properties  of  the  solar  cell  affect  this  energy  conversion? 

3.  What  are  the  advantages  and  disadvantages  of  concentrated  solar  power 
systems  over  conventional  photovoltaic  systems? 

4.  At  what  level  do  the  cunent  physical  models  predict  the  behavior  of 
concentrated  solar  power  systems? 

5.  How  can  the  existing  methods  be  improved  to  develop  a  new  dynamic  thermal 
model  for  a  concentrated  solar  power  system? 

6.  What  benefits  do  the  model  presented  provide  to  the  existing  technology? 

7.  How  could  this  research  be  expanded  in  the  future? 

8.  What  are  the  strategic  implications  of  renewable  energy  use  by  the  military? 
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The  increased  understanding  of  concentrated  PV  systems  gleaned  from  this 
research  will  be  very  beneficial  to  the  U.S.  Air  Force  and  Department  of  Defense.  As 
with  most  civilian  technologies,  further  military  research  into  terrestrial  PV  systems  will 
most  certainly  drive  the  development  of  these  technologies  ahead  farther  and  faster  in  the 
future. 
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II:  Theory 


The  theoretical  analysis  of  a  semiconductor  solar  cell  must  begin  with  a 
fundamental  understanding  of  the  solid  state  physics  governing  this  device.  The  first  goal 
of  this  analysis  is  to  investigate  the  simple  principle  of  a  PN -junction  under  thermal 
equilibrium,  which  operationally  represents  a  solar  cell  in  the  dark.  Several  important 
tenns  are  presented  and  defined  in  order  to  lay  the  framework  for  the  following  analysis. 
The  important  dependency  upon  the  excess  minority  carrier  concentration  in  the  device  is 
described.  Equilibrium  diode  currents  and  their  associated  equations  are  also  introduced 
and  applied  to  the  analysis.  The  optical  photon  absorption  of  semiconductors  is 
investigated  introducing  the  photogenerated  current  produced  across  the  PN -junction. 
Next,  the  nonequilbrium  process  of  carrier  generation  caused  by  solar  illumination  is 
discussed,  as  applied  to  the  simple  PN-junction.  Building  upon  this  analysis,  the  non¬ 
uniformity  of  carrier  generation  is  discussed  along  with  the  role  that  recombination  plays 
in  each  region  of  the  cell.  Synthesizing  the  analysis,  an  expression  for  the  excess 
minority  carrier  concentrations  on  each  side  of  the  junction  and  the  resulting  total  current 
through  the  solar  cell  is  derived.  Finally,  the  fundamental  solar  cell  parameters  of  open 
circuit  voltage  Voc,  short  circuit  current  Isc,  and  fill  factor^  are  defined  in  the  context  of 
the  preceding  analysis. 

2,1  Semiconductor  Physics  Background 

The  physical  principles  describing  semiconductor  solar  cells  are  rooted  in 
quantum  theory,  the  basics  of  which  explain  the  concept  of  a  semiconducting  material.  A 
semiconductor  is  defined  by  the  Science  Dictionary  as  a  material  with  conductive 
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properties  between  those  of  a  conductor  and  an  insulator,  and  semiconductors  are 
characterized  by  their  atomic  crystal  structure  [10],  Figure  3  shows  a  plot  of  several 
materials  based  on  their  electrical  conductivities  [8].  Referring  to  the  figure,  common 
insulators  are  on  the  left,  semiconductors  are  in  the  middle,  and  conductors  are  on  the 
right. 
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Figure  3:  Semiconductor  Electrical  Properties  [8] 

Referring  to  Figure  3  another  fundamental  property  of  semiconductors  is 
revealed-  the  variable  electrical  conductivity  (and  thus  resistivity).  Therefore,  the 
electrical  properties  of  semiconductors  can  be  effectively  “tuned”  to  generate  the  specific 
values  of  impedance  over  a  given  range  depending  on  the  temperature  and  especially  the 
impurity  doping  properties  of  the  device. 

The  crystal  structure,  or  lattice,  of  a  semiconductor  is  extremely  important  to  the 
electrical  behavior  of  the  material.  All  semiconducting  materials  are  found  on  the  right 
side  of  the  periodic  table. 


10 


VIIIA 

2 


29 


IB  MB 

130.  ‘  ■ 

Zn 


Cu 

63.54 


65.37 


Ag  Cd 

107.870  112.40 

79  80 

Au 

196.967  200. 


III  A 

IVA 

VA 

VIA 

VI 1 A 

He 

4.003 

5 

B 

6 

C 

N 

8 

0 

9 

F 

10 

Ne 

10.811 

12.011 

14.007 

15.999 

18.998 

20.183 

13 

Al 

14 

Si 

15 

P 

16 

s 

17 

Cl 

18 

Ar 

26.982 

28.086 

30.974 

32.064 

35.453 

39.948 

31 

Ga 

32 

Ge 

33 

As 

Se 

35 

Br 

36 

Kr 

69.72 

| 72.59 

74.922 

78.96 

79.909 

83.80 

49 

In 

50 

Sn 

Sb 

52 

Te 

53 

1 

54 

Xe 

114.82 

118.69 

121.75 

127.60 

126.904 

131.30 

81 

Tl 

82 

Pb 

83 

Bi 

84 

Po 

85 

At 

86 

Rn 

204.37 

207.19 

208.980 

(210) 

(210) 

(222) 

Figure  4:  Right  Side  of  the  Periodic  Table  Showing  the  Semiconducting  Elements 

Most  common  semiconducting  materials  are  found  in  groups  III,  IV,  and  V, 
which  is  very  significant  because  the  periodic  group  of  elements  is  organized  by  outer 
valence  electron  structure,  where  group  III  elements  have  three  valence  electrons, 
elements  in  group  IV  have  four  valence  electrons,  and  accordingly  elements  in  group  V 
have  five  valence  electrons. 

Semiconductors  are  found  in  both  elemental  (e.g.  Si)  and  compound  fonn  (e.g. 
GaAs),  but  every  semiconductor  material  is  characterized  by  the  properties  of  its  crystal 
lattice.  The  covalent  bonds  formed  within  a  semiconducting  material  determine  the 
shape  of  the  crystal  lattice  [8].  For  an  in  depth  explanation  on  the  quantum  physics 
governing  this  interaction,  a  Solid  State  Physics  text  book  such  as  McKelvey’s  is  desired, 
but  for  the  sake  of  this  analysis  applied  to  solar  cells,  a  cursory  discussion  is  sufficient 
[10].  Most  semiconductors  are  composed  of  the  diamond  lattice,  where  each  atom  has  a 
covalent  bond  with  four  adjacent  atoms.  Compound  semiconductors  such  as  GaAs  have 
a  crystal  lattice  similar  to  the  diamond  lattice,  but  since  the  crystal  is  comprised  of  two 
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different  elements,  the  compound  semiconductor  crystal  lattice  is  referred  to  as  the  zinc 
blende  structure  as  shown  in  Figure  5. 
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Figure  5:  Diamond  and  Zinc  Blende  Lattice  Structure  [10] 

With  an  appropriate  understanding  of  the  semiconductor  crystal  structure  in  hand, 
the  electrical  properties  of  these  materials  can  be  investigated. 

2.2  Electrical  Properties  of  Semiconductors 

In  crystalline  materials  such  as  semiconductors  the  electrical  properties  are  driven 
by  different  energy  levels.  Electrons  within  a  semiconductor  are  restricted  to  certain 
energy  levels  (bands)  between  the  ground  state,  in  which  the  electrons  are  tightly  bound 
to  the  nucleus  of  each  atom,  and  the  free  state,  in  which  electrons  are  readily  available  for 
conduction  [8].  The  two  most  important  energy  bands  within  a  semiconductor  material 
are  the  valence  band  and  the  conduction  band.  The  valence  band  Ev  is  energetically  close 
to  the  nucleus  and  nearly  full  of  electrons,  while  the  conduction  band  Ec  (next  energy 
level  up  from  the  valence  band)  is  only  slightly  filled  with  electrons.  The  amount  of 
energy  necessary  to  excite  an  electron  from  the  valence  band  to  the  conduction  band  is 
the  band  gap  energy  Eg,  and  physically  describes  the  energy  difference  between  these  two 
bands  [9].  Semiconductor  materials  can  be  generally  characterized  as  direct  band  gap 
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materials  or  indirect  band  gap  materials.  This  characterization  is  based  upon  the  nature 
of  the  energy  levels  within  the  material.  Both  the  valence  and  conduction  bands  can  be 
described  as  wave-like  energy  levels  with  varying  momentum.  In  direct  band  gap 
materials,  the  band  gap  energy  corresponds  to  a  wave  k-vector  equal  to  zero  where  the 
minimum  distance  between  the  valence  and  conduction  bands  occurs  at  the  same  value  of 
electron  momentum,  as  shown  in  Figure  6  [10]. 


Figure  6:  Direct  Band  Gap  Semiconductor  [10] 

In  indirect  band  gap  materials,  the  k-vector  does  not  equal  zero  and  the  minimum 
distance  between  the  two  bands  occurs  indirectly  and  corresponds  to  a  difference  in 
electron  momentum  between  the  two  bands.  Therefore,  in  an  indirect  band  gap  material 
shown  in  Figure  7,  an  electron  cannot  travel  to  the  conduction  band  without  a 
corresponding  change  in  momentum  which  requires  a  phonon  assisted  transition,  also 
known  as  a  lattice  vibration  [10]. 
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Figure  7:  Indirect  Band  Gap  Semiconductor 


The  intermediate  energy  level  located  between  the  valence  and  conduction  bands 
is  known  as  the  Fermi  level  Ef,  and  is  often  used  as  a  reference  level  in  the  analysis  of 
semiconductors  as  will  be  shown  in  the  following  analysis.  The  diagram  in  Figure  8 
shows  the  standard  representation  of  an  intrinsic  (pure)  semiconductor  material,  including 
the  three  energy  levels  discussed  above. 


Figure  8:  Equilibrium  Energy  Band  Diagram 
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It  is  standard  and  accepted  practice  to  draw  the  band  diagram  of  direct  and 
indirect  band  gap  materials  identically,  since  the  physical  significance  of  this  difference 
has  primarily  one  consequence  to  be  discussed  later  in  the  analysis. 

Classically,  there  are  two  main  charge  carriers  within  a  semiconductor:  electrons 
and  holes,  and  each  of  these  carriers  contribute  to  the  total  electricity  generated  by  a 
photovoltaic  cell.  A  semiconductor  material  with  an  excess  concentration  of  electrons  is 
referred  to  as  an  //-type  semiconductor,  whereas  a  semiconductor  with  an  excess 
concentration  of  holes  is  referred  to  as  p- type.  Furthermore,  in  extrinsic  semiconductors 
the  Fenni  level,  which  physically  describes  the  chemical  potential  of  the  material,  is 
closer  to  the  conduction  band  Ec  for  //-type  materials  and  closer  to  the  valence  band  Ev  for 
/j-type  materials. 

2.3  Equilibrium  P-N  Junction 

The  elemental  device  structure  of  most  modern  solar  cells  is  the  PN-junction,  in 
which  a  semiconductor  material  of  one  conductivity  type  (n  or  p)  is  electrostatically 
joined  with  a  semiconductor  material  of  another  type  (n  or  p).  The  conductivity  type  of 
the  semiconductor  material  is  driven  by  the  impurity  doping  concentration  on  each  side 
of  the  junction,  where  //-type  semiconductors  have  an  excess  of  electrons  and  p- type 
semiconductors  have  an  excess  of  holes.  Extrinsic  doping  changes  the  conductivity  of 
the  semiconductor  by  adding  charge  carriers,  resulting  in  a  change  in  Ef,  either  raising  or 
lowering  the  Fermi  energy  of  the  semiconductor  relative  to  its  intrinsic  value.  A  PN- 
junction  under  thermal  equilibrium  is  shown  in  Figure  9. 
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Figure  9:  Equilibrium  PN-junction  Band  Diagram 

When  the  two  sides  of  the  PN-junction  are  joined  together,  the  difference  in  the 
intrinsic  Fenni  energies  creates  a  built  in  potential  Vy,  also  known  as  the  potential 
barrier,  across  the  junction  but  the  total  Fermi  energy  level  across  the  entire  device 
remains  constant  as  shown  in  Figure  9.  In  order  to  conduct  charge  across  the  junction, 
the  barrier  height  must  be  exceeded  by  an  externally  applied  voltage  Va.  Furthermore,  at 
the  junction  a  depletion  region  (also  known  as  a  space  charge  region)  forms  on  either  side 
where  no  carriers  are  present.  The  depletion  region  arises  because  of  carrier  diffusion  at 
the  junction,  and  the  electric  field  over  this  region  causes  these  carriers  to  move  from  one 
side  of  the  depletion  region  to  the  other,  provided  the  carriers  are  located  within  a 
diffusion  length  of  the  depletion  region.  The  importance  of  this  phenomenon  cannot  be 
understated,  since  it  is  the  electric  field  over  the  depletion  region  that  drives  the  main 
process  by  which  light  is  converted  into  electricity.  However,  analytical  investigation  of 
this  electric  field  involves  a  very  important,  fundamental  assumption. 
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2.4  Poisson’s  Equation 

The  energy  band  diagram  gives  only  a  qualitative  description  of  how  the  electric 
field  varies  over  the  junction,  but  Poisson’s  equation  can  be  used  to  analytically  describe 
the  charge  density,  electric  field,  potential,  and  depletion  region  width  of  a  PN-junction. 
Poisson’s  equation  is  written  below,  where  (p  is  the  potential,  p  is  the  charge  density,  q  is 
the  elemental  charge,  n  and  p  are  the  electron  and  hole  densities,  Nd+  and  Na"  are  the 
ionized  electron  and  hole  densities,  and  ssc  is  the  relative  pennittivity  of  the 
semiconductor  [14]. 


d2w  p  q 

-^=~r=~  J-(P-n  +  Na-Na) 

L/A,  C.sc  Csc 


(2-1) 


The  depletion  approximation  allows  the  complex  equations  for  the  electric  field, 
potential,  and  depletion  region  widths  to  be  numerically  solved.  The  depletion 
approximation  assumes  that  no  carriers  are  present  in  the  depletion  region,  thus  this 
region  is  fully  depleted  of  charge  carriers.  The  depletion  approximation  also  assumes 
that  the  depletion  region  has  well  defined  boundaries,  so  that  the  transition  from  the 
depletion  region  to  the  bulk  n  and  p  (or  quasi  neutral  regions)  regions  is  abrupt.  The 
justification  for  use  of  the  depletion  approximation  in  the  electrostatic  analysis  of  a  PN- 
junction  stems  from  the  fact  that  carrier  concentrations  change  exponentially  with  the 
Fermi  energy  relative  to  the  conduction  band  edge,  and  this  change  occurs  over  the 
depletion  region  as  shown  in  Figure  9  [12].  Accordingly,  the  contributions  of  n  and  p 
over  the  depletion  region  become  negligible  and  can  be  removed  from  the  above 
equation,  simplifying  the  following  analysis. 
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dx2=-^c  =  -^-^ 


(2-2) 


The  charge  density  can  then  be  derived  for  four  distinct  regions  of  the  device. 

p(x)  =  0  xn  <  x  <  oo  (2-3.1) 

p(x)  =  —  qNa  0  <  x  <  xn  (2-3.2) 

p(x)  =  qNd  —xp  <  x  <  0  (2-3.3) 

p(x)  =  0  —  oo  <  x  <  — xp  (2-3.4) 

A  plot  of  the  charge  density  as  a  function  of  doping  and  distance  in  a  PN-junction 
is  shown  below  in  Figure  10. 


Figure  10:  PN-junction  Charge  Density  Profile  [13] 

— 1  r/  F  (  y) 

Since  charge  is  related  to  electric  field  E  by  ^  =  q(Nd  —  Na ),  the  expressions 

for  charge  on  each  side  of  the  junction  can  be  integrated  to  yield  the  electric  field  over  the 

depletion  region  as  a  function  of  doping  and  distance  into  the  material. 

qNa(x  +  x„ 

£(X)  =  - E 

£qr 


0  <  X  < 


(2-4.1) 
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£(x)  = 


(2-4.2) 


qNd(x-xn) 


—xp  <  x  <  0 


A  plot  of  the  electric  field  over  the  depletion  region  reveals  that  the  greatest 
magnitude  of  the  field  occurs  at  the  PN-junction  as  shown  in  Figure  11. 


Figure  11:  Equilibrium  PN-junction  Electric  Field  [13] 


Again,  the  importance  of  this  electric  field  over  the  depletion  region  cannot  be 
understated.  As  will  be  shown  in  the  following  analysis,  the  electron  drift  current  over 
this  region  complements  the  diffusion  current  to  drive  the  electrical  output  of  PN-junction 
solar  cells.  Further  iterations  of  Poisson’s  equation  can  be  solved  to  derive  expressions 
for  built  in  potential  cp;  (Vbi)  and  depletion  region  width  W  in  the  device,  where  Va  is  the 
applied  voltage  (zero  under  thermal  equilibrium).  The  result  of  these  calculations  are 
shown  in  equations  (2-5)  and  (2-6),  but  an  in  depth  description  of  each  step  is  given  by 
Muller  and  Kamins  [14], 


(2-5) 
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w  = 


(2-6) 


With  expressions  for  the  charge,  electric  field,  potential,  and  depletion  region 
width  derived,  the  physical  movement  of  charge  carriers  within  an  equilibrium  PN- 
junction  must  be  investigated. 

2.5  Equilibrium  Charge 

Two  equal  and  opposite  forces  dominate  the  movement  of  charges  within  an 
equilibrium  PN -junction.  Electron  drift  is  the  first  mechanism  that  causes  carriers  to 
move  in  a  PN-junction  and  is  caused  by  the  electric  field  created  over  the  depletion 
region.  The  electric  field  causes  carriers  with  the  opposite  charge  to  “drift”  to  each  side 
of  the  junction  and  with  respect  to  the  band  diagram  shown  below  in  Figure  12  causes 
electrons  to  roll  downhill  from  the  p-side  to  the  n-side  and  holes  to  go  uphill  from  the 
opposite  direction. 
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Figure  12:  Charge  Movement  Within  Equilibrium  PN-junction  [16] 
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Under  thennal  equilibrium,  the  process  of  electron  drift  is  exactly  balanced  by  the 
process  of  electron  diffusion  which  is  driven  by  the  concentration  gradient  of  holes  and 
electrons  at  the  PN-junction  [11].  In  much  the  same  way  that  a  gas  diffuses  into  air  and 
spreads  out  into  a  given  volume,  electrons  and  holes  diffuse  throughout  the  device  once 
the  two  regions  are  brought  into  contact.  Additionally,  since  the  minority  carrier 
concentrations  on  the  n  and  /;-sides  will  be  greatest  at  the  edge  of  the  depletion  regions, 
the  concentration  gradient  causes  minority  carriers  to  diffuse  away  from  the  depletion 
boundaries  towards  the  edges  of  the  semiconductor  material.  The  distance  that  the 
minority  carriers  diffuse  into  the  semiconductor  material  is  defined  by  a  parameter  known 
as  the  diffusion  length  Ln/P .  As  will  be  shown  in  the  following  analysis,  once  the  PN- 
junction  is  disturbed  from  equilibrium  electron  diffusion  will  be  the  primary  mechanism 
that  drives  the  movement  of  charge  within  the  bulk  of  the  semiconductor  material. 

As  stated  in  the  introduction  of  this  analysis,  the  equilibrium  PN-junction 
operationally  represents  a  solar  cell  in  the  dark-  no  net  current  flows,  no  external  bias  is 
applied,  and  the  excess  minority  carrier  concentrations  are  negligible.  However,  to  make 
this  analysis  useful  the  condition  of  solar  illumination  must  be  applied,  thus  disturbing 
the  simple  PN-junction  from  equilibrium  and  making  the  analysis  much  more  interesting. 

2.6  Solar  Illumination 

Before  delving  into  the  details  of  the  nonequilbrium  PN-junction  analysis  it  is 
beneficial  to  explain  some  of  the  important  fundamental  properties  of  solar  illumination. 
A  solar  cell  converts  sunlight  into  electricity  by  means  of  the  photovoltaic  effect  [14]. 
When  light  falls  upon  the  surface  of  a  semiconductor,  incident  photons  with  enough 
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energy  break  the  covalent  bonds  within  the  material  and  free  up  the  valence  electrons  for 
conduction.  The  photovoltaic  effect  is  rooted  in  a  concept  known  as  the  wave-particle 
duality  of  light,  which  states  that  light  can  simultaneously  be  treated  as  an 
electromagnetically  radiated  wave  and  as  a  particle  with  finite  energy  [8].  The  latter 
definition  lends  itself  well  to  the  solid  state  analysis  of  solar  cells. 

The  energy  in  solar  radiation  is  given  by  E  —  hv  where  h  is  Plank’s  constant  and 
v  is  the  frequency  of  the  light.  Figure  13  shows  AM  1.5  solar  irradiance  spectrum  that  is 
used  as  the  standard  for  analysis  in  most  terrestrial  solar  energy  systems.  The  frequencies 
of  light  most  irradiant  at  the  surface  of  the  earth  range  from  300-750  THz.  Accordingly 

C 

since  frequency  is  related  to  wavelength  by  1  =  -,  where  c  is  the  constant  speed  of  light 

o 

(3x10  m/s),  a  common  form  for  the  wavelength  dependent  photon  energy  of  light 

1  24 

is  E  —  — — ,  where  X  is  in  pm  [14].  The  wavelengths  of  light  falling  on  the  cell  above  the 
band  gap  energy  will  vary  from  .4-1  pm  and  the  associated  energy  range  is  1. 2-3.1  eV. 
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TYPICAL  SOLAR  SPECTRUM 


Figure  13:  AM  1.5  Solar  Spectrum 

The  bandgap  energy  is  the  amount  of  energy  required  to  free  an  electron  from  its 
covalent  bond  in  the  semiconductor  material,  thus  allowing  it  to  contribute  to  the  transfer 
of  charge  (electricity).  The  photon  energy  of  the  light  must  be  at  least  as  great  as  the 
bandgap  energy  of  the  semiconductor  in  order  to  free  an  electron  and  generate  a  current 
[15].  As  incident  photons  with  energy  greater  than  (or  equal  to)  the  bandgap  energy  of 
the  semiconductor  solar  cell  fall  onto  the  PN-junction,  electrons  are  excited  and  electron- 
hole  pairs  are  generated.  However,  photon  energies  greater  than  the  bandgap  energy  also 
generate  heat  leading  to  thermal  losses  as  the  excited  electron  falls  back  to  the  edge  of  the 
conduction  band.  The  goal  of  the  solar  cell  designer  is  to  closely  match  the  bandgap 
energy  of  the  semiconductor  cell  with  the  photon  energy  of  incident  light  so  that  the 
energy  production  can  be  maximized  while  minimizing  losses  due  to  heat. 

2.7  Photon  Absorption  and  Generation 
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The  most  fundamental  mechanism  in  the  photovoltaic  energy  conversion  process 
is  photon  absorption  and  the  associated  electron-hole  pair  (EHP)  generation.  Most  solar 
cell  models  assume  that  every  photon  absorbed  by  the  semiconductor  material  generates 
exactly  one  EHP  [9].  Previously,  it  was  explained  that  photons  with  energy  below  the 
semiconductor  bandgap  energy  would  not  contribute  to  conduction  and  were  assumed  to 
transmit  entirely  through  the  material.  However,  the  absorption  of  photons  is  a  statistical 
process,  and  it  must  be  stated  that  some  photons  below  the  band  gap  will  contribute  to  the 
total  heat  generated  within  the  solar  cell,  without  directly  contributing  to  the  conduction 
of  electricity  across  the  PN-junction.  Furthermore,  not  all  photons  with  energy  above  the 
band  gap  will  contribute  to  conduction.  The  quantum  efficiency  [3  describes  the 
percentage  of  photons  incident  on  to  the  cell  that  actually  contribute  to  charge 
conduction.  The  details  of  this  phenomenon  will  be  discussed  later  in  the  analysis  and 
applied  to  the  dynamic  thermal  model.  Regardless,  for  the  electrical  model  of  the  PN- 
junction  solar  cell,  the  former  assumption  that  every  photon  absorbed  generates  one  EHP 
will  be  used. 

It  is  important  to  realize  that  photon  absorption  and  electron  generation  are 
wavelength  dependent  properties,  since  different  wavelengths  of  light  carry  different 
amounts  of  photon  energy.  Shorter  wavelengths  of  light  with  higher  photon  energies, 
blue  light  for  example,  will  be  absorbed  more  readily  by  the  material  and  thus  do  not 
travel  as  far  into  the  PN-junction  solar  cell  [17].  It  can  be  assumed  that  for  shorter 
wavelengths,  almost  all  of  the  generation  occurs  near  the  surface  of  the  cell.  However,  at 
longer  wavelengths  of  light  with  lower  photon  energies,  red  light  for  example,  the 
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photons  penetrate  on  average  entirely  though  the  bulk  of  the  device.  Furthermore,  due  to 
the  physical  structure  of  the  solar  cell,  all  of  the  carrier  generation  will  occur  over  the  two 


diffusion  lengths  on  the  n  and  p-side  and  the  entire  depletion  region,  or  Ln  +  W  +  lp.  A 


pictorial  interpretation  of  this  is  shown  in  Figure  14. 
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Figure  14:  Photon  Absorption  at  the  PN-junction  [15] 


The  wavelength  (and  thus  frequency)  dependent  photon  flux  is  given  by  S(x,  A)  — 
S0(A)e~ax  where  So  is  the  wavelength  dependent  number  of  incident  photons  and  a  is  the 
absorption  coefficient,  a  material  dependent  parameter  that  describes  how  a  given 
material  generates  electrons  in  a  device.  Figure  15  shows  the  absorption  coefficient  as  a 
function  of  wavelength  at  300K  for  7  different  semiconductor  materials: 
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Figure  15:  Wavelength  Dependent  Photon  Absorption  for  Different  Semiconductor  Materials  [16] 


Since  generation  is  dependent  upon  photon  absorption,  the  generation  rate  is  also 
a  function  of  wavelength.  However,  it  is  useful  to  sum  the  wavelength  dependent 
magnitude  of  solar  irradiation  in  order  to  investigate  how  the  generation  rate  varies  as  a 
function  of  distance  into  the  material  over  the  entire  spectrum.  Figure  16  shows  a 
diagram  of  the  net  generation  rate  versus  position  for  silicon. 
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Figure  16:  Generation  Rate  for  Si 

Another  important  optical  property  to  discuss  in  the  photovoltaic  conversion  of 

energy  is  the  reflection  of  light  off  of  the  semiconductor  surface.  Often  times,  the 

reflectivity  R  is  a  significant  detractor  to  the  overall  photon  absorption  so  this 

phenomenon  cannot  be  ignored  when  analyzing  the  generation  rate  within  a 

semiconductor.  The  fraction  of  incident  light  that  is  reflected  is  a  function  of  the  index  of 

refraction  n  and  the  extinction  coefficient  k  of  the  material  [16]. 

(n  -  l)2  +  k2  (2-7) 

(n  +  l)2  +  k2 

With  photon  absorption  and  reflection  accounted  for,  a  more  realistic  wavelength- 
dependent  description  of  carrier  generation  is  revealed,  where  G(x,  /.)  denotes  position 
and  wavelength  dependent  generation. 

G{x,X )  =  /?(A)a(A)S0(A)(l  -  R)e~a^x  (2-8) 
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For  most  homojunction  solar  cells,  the  quantum  efficiency  P  approaches  unity  for 
all  wavelengths  with  energies  greater  than  the  bandgap  energy,  and  for  this  analysis  the 
net  generation  is  rewritten  with  Go  as  the  initial  generation  at  the  surface  accounting  for 
the  reflectivity  and  number  of  incident  photons. 

G  =  60e(_ax)  (2-9) 

The  important  points  to  remember  are  that  generation  is  greatest  at  the  surface  and 
decreases  exponentially  throughout  the  bulk  of  the  cell.  As  was  shown  previously,  the 
absorption  coefficient  is  wavelength  dependent  so  it  follows  that  generation  is  a 
wavelength  dependent  process  as  well. 

2.8  Carrier  Recombination 

An  equally  important  process  to  carrier  generation  that  must  be  discussed  is 
carrier  recombination.  Recombination  serves  to  counteract  generation,  in  order  to  return 
the  nonequilibrium  cell  to  equilibrium,  and  physically  describes  the  recombining  of  holes 
and  electrons  in  the  cell.  Once  an  electron  and  hole  recombine,  the  neutral  charge  is  no 
longer  available  for  conduction;  thus,  it  is  important  to  minimize  carrier  recombination 
throughout  the  cell  in  order  to  maximize  generation  and  charge  conduction.  The  process 
of  recombination  is  most  prevalent  at  impurities  or  defects  within  the  crystal  structure  of 
the  cell  and  also  at  the  surface  of  the  cell.  Accordingly,  since  the  process  of  doping 
intentionally  introduces  impurities  to  the  crystal  structure  to  change  the  electrical 
properties  of  the  cell,  doping  consequently  increases  recombination  in  the  cell  [16]. 
Furthennore,  immaculate  semiconductor  processing  is  paramount  since  any  unintentional 
impurities  will  only  enhance  the  carrier  recombination. 
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There  are  three  types  of  recombination  in  the  bulk  of  a  semiconductor  solar  cell: 
radiative  recombination,  Auger  recombination,  and  Shockley-Read-Hall  (SRH) 
recombination.  Radiative  recombination,  also  referred  to  as  band-to-band  recombination, 
describes  the  phenomenon  where  an  electron  directly  recombines  with  a  hole  in  the 
conduction  band  releasing  a  photon  of  light.  The  process  of  radiative  recombination 
drives  the  operation  of  light  emitting  diodes  (LED)  and  semiconductor  lasers  but  is 
detrimental  to  the  operation  of  a  photovoltaic  cell.  Radiative  recombination  is  usually 
neglected  in  the  analysis  of  indirect  bandgap  semiconductors  such  as  Silicon  (Si),  but  this 
process  is  potentially  significant  for  direct  bandgap  semiconductors  such  as  Gallium 
Arsenide  (GaAs). 

Auger  recombination  is  another  unavoidable  recombination  mechanism  that 
describes  the  phenomenon  where  an  electron  and  hole  recombine,  and  give  the  resulting 
excess  kinetic  energy  to  another  electron  in  the  conduction  band.  The  excess  energy 
given  to  the  third  carrier  in  the  conduction  band  is  expended  in  the  cell  as  heat.  Auger 
recombination  dominates  in  concentrated  solar  cells  [21]. 

The  final  recombination  mechanism  in  semiconductor  solar  cells  is  SRH 
recombination.  SRH  recombination  is  due  to  crystal  defects  and  is  thus  most  important 
from  a  processing  perspective,  and  is  usually  the  most  prevalent  recombination  process  in 
a  solar  cell.  Crystal  defects,  whether  implanted  intentionally  or  the  product  of  imperfect 
processing  techniques,  introduce  energy  states  within  the  bandgap  energy  of  a  PN- 
junction.  These  pseudo-energy  levels  known  as  “traps”  are  recombination  centers  in  the 
bandgap  since  any  electron  that  moves  from  the  valence  band  to  a  trap  is  likely  lost  to 
recombination,  especially  if  the  trap  is  near  the  mid-band  energy  of  the  semiconductor. 


29 


SHR  recombination  is  a  two-step  process  where  the  electron  moves  from  the  valence 
band  to  the  trap,  and  then  either  up  to  the  conduction  band  or  back  to  the  valence  band 
where  it  recombines.  The  process  is  functionally  the  same  but  reversed  for  holes.  The 
rate  of  SHR  recombination  depends  on  the  trap  energy  levels,  where  traps  near  the  band 
edges  are  less  likely  to  lead  to  recombination. 

Two  material  parameters  that  are  affected  by  recombination  processes  are  carrier 
lifetime  x  and  diffusion  length  D.  The  carrier  lifetime  physically  describes  the  amount  of 
time  an  EHP  exists  in  a  PN-junction  after  generation  and  before  recombination  (or 
conduction).  Therefore,  longer  carrier  lifetimes  are  associated  with  less  carrier 
recombination.  The  diffusion  length  describes  the  physical  distance  an  electron  or  hole 
travels  before  recombining.  Ideally,  the  diffusion  length  is  maximized  to  reduce 
recombination.  The  importance  of  these  parameters  with  respect  to  the  electrical  output 
of  a  PN-junction  solar  cell  will  be  discussed  later  in  the  analysis.  Furthermore,  each  of 
the  three  recombination  processes  will  be  treated  simultaneously  but  it  is  beneficial  to 
understand  the  three  different  mechanisms  that  an  electron  and  hole  recombine  in  order  to 
minimize  the  detrimental  effect  to  the  overall  efficiency  and  electrical  output. 

2.9  PN-junction  Solar  Cell  Current 

With  an  understanding  of  the  equilibrium  PN-junction  coupled  with  an 
understanding  of  the  role  of  light  in  photovoltaic  energy  conversion,  it  is  necessary  to 
investigate  how  the  PN-junction  changes  under  nonequilibrium  conditions.  To  facilitate 
this  analysis  the  ideal  diode  equation  is  applied  to  the  solar  PN-junction  and  the  total 
current  due  to  solar  illumination  is  investigated.  The  value  of  the  following  model  is  that 
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the  exponential  generation  of  charge  carriers  is  incorporated  into  the  diode  current 
derived  from  the  transport  and  continuity  equations.  Most  every  solar  cell  current  model 
currently  in  the  literature  refers  to  the  familiar  current  equation  that  breaks  the  diode 
current  ID  and  photogenerated  current  Iph  into  two  separate  parts,  where  Itot  is  the  total 
current  flowing  through  the  cell: 


hot  —  Ipi i 


(2-10) 


The  photogenerated  current  is  commonly  given  by  Iph  —  qG0A{Ln  +  W  +  Lp), 

1  —  1 

and  the  diode  current  is  given  by  Id  —  I0  e^r  —  1 1,  where  I0  is  found  using  the  transport 
and  continuity  equations  [12]. 


2.10  Solar  Cell  Current-Voltage  Relationship 


A  useful  method  to  qualitatively  examine  the  operation  of  a  solar  cell  over  the 
entire  range  of  operating  conditions  is  to  plot  the  solar  cell  current  versus  the  voltage, 
known  as  the  IV  curve.  An  example  IV  curve  is  shown  in  Figure  17. 
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Solar  Cell  IV-tuive 


Figure  17:  Solar  Cell  IV  Curve 

The  IV  characteristics  of  the  cell  clearly  show  the  short  circuit  current  Isc  and  the 
open  circuit  voltage  Voc,  the  y-intercept  and  x-intercept  respectively,  and  reveal  another 
important  solar  parameter  known  as  the  fill  factor,  denoted  by  ff.  The  fill  factor  describes 
the  “squareness”  of  the  IV  curve  since  the  ideal  current-voltage  relationship  is 
rectangular.  Therefore,  it  is  beneficial  to  have  the  fill  factor  as  close  as  possible  to  unity 
but  this  can  never  be  achieved  due  to  the  exponential  nature  of  PN -junction  properties. 

An  expression  for  the  fill  factor  is  defined  below,  where  Im  is  the  maximum  operating 
current  and  Vm  is  the  maximum  operating  voltage. 

IjrXn  (2-11) 

n  IscVoc 

The  maximum  operating  conditions  Im  and  Vm  are  the  values  where  rectangular 
area  within  the  IV-curve  is  maximized.  Furthermore,  the  magnitude  of  the  fill  factor 
depends  upon  the  device  structure  and  material  properties,  but  most  fill  factor  values 
typically  range  from  0.7-0.85  [14]. 
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2.11  Load  Resistance 


The  explanation  of  fill  factor  leads  into  another  important  component  of  any  solar 
power  system,  the  load.  It  is  necessary  to  connect  the  solar  cell  to  a  load  in  order  to 
generate  power.  While  the  shape  of  the  IV  curve  is  independent  of  the  load,  the  operating 
point  (Iop,  Vop)  is  determined  from  the  load  resistance.  Accordingly,  the  operating  point 
can  be  tuned  to  get  the  current  and  voltage  as  close  as  possible  to  Im  and  Vm  by  matching 
the  load  resistance  to  the  cell.  Load  (impedance)  matching  is  an  extremely  important 
concept  in  the  physical  set  up  of  solar  arrays. 

The  operating  point  of  the  solar  cell  based  on  given  load  resistance  Rl  is  most 

1 

easily  found  using  a  load  line  construction,  a  linear  line  with  slope  — .  Figure  18  below 

rl 

shows  a  load  line  using  a  1 .68  Q  load,  for  the  solar  cell  used  in  the  previous  IV  curve 
figure. 


Solar  Cell  IV-curve 


Figure  18:  Load  Line  Plotted  on  IV-Curve 
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When  the  cell  is  connected  to  a  load  the  voltage  across  and  current  through  the 
load  are  equal  to  the  voltage  and  current  produced  by  the  cell,  but  the  current  flows  in  the 
opposite  direction.  Using  Ohm’s  law,  an  expression  for  the  current  through  the  load  can 


be  written. 


fop  — 


Kp 

Rl 


(2-12) 


2,12  Power  and  Efficiency 


The  load  serves  to  transform  the  electrical  output  of  the  cell  into  usable  energy  in 
the  fonn  of  power.  Power  is  the  fundamental  unit  of  electrical  output  for  most  energy 
systems,  and  most  every  solar  array  is  described  in  tenns  of  the  total  power  generated. 
Again,  this  power  is  dependent  upon  the  load  and  is  optimized  under  load-matched 
conditions.  The  relationship  between  voltage,  current  and  power  is  fundamental  in 
electronics. 


P  =  IV  (2-13) 

The  most  practical  measure  of  merit  for  solar  cells  is  the  electrical  efficiency  q,  a 
ratio  of  the  power  incident  upon  the  cell  in  the  form  of  light  Pm  to  the  power  generated  by 
the  cell  Pout. 


V  = 


out 


x  100 


(2-14) 


Solar  cell  efficiencies  are  generally  in  the  range  of  12%  for  most  terrestrial  non¬ 
concentrated  solar  cells  up  to  25%  for  highly  specialized  multi -junction  cells  designed  for 
use  in  space-based  or  CSP  systems.  Several  factors  contribute  to  efficiency  losses  in 
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solar  power  systems,  and  Figure  19  summarizes  the  major  factors  decreasing  this 
electrical  efficiency. 

1 00%  Incident  radiation 

Insufficient  photon  energy 
hv  <Eg 

Excessive  photon  energy 
Near  suffice  EHP  recombination 
hv  >  Eg 

Collection  efficiency  ofphotons 
Voc  *  (0.6 Egy(ekB) 

FF  *  0.85 


A  noAiintirur  for  i;ori oi i c  Ioocpc  r\f  oriorm;  in  a  Vii  rrlo  offioionov  Qi 

Figure  19:  Detractors  to  Overall  Efficiency  [16] 

With  an  understanding  of  the  major  solar  cell  measures  of  merit,  the  theoretical 
analysis  of  the  simple  PN-junction  solar  cell  is  complete. 

2.13  Theoretical  Synopsis 

To  recap  the  analysis,  the  theoretical  discussion  of  the  PN-junction  solar  cell 
began  with  an  investigation  of  the  material  properties  of  semiconductors  and  an 
explanation  of  their  characteristic  crystal  structure.  A  short  explanation  of  the  atomic 
bonding  within  a  semiconductor  material  was  discussed  and  the  periodic  table  was 
referenced  showing  groups  III-V,  the  most  important  groups  to  the  semiconductor 
industry. 

Following  the  atomic  level  discussion,  an  introduction  to  the  quantum  physics 
governing  the  electrical  properties  of  semiconductors  was  presented.  Several 
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fundamental  properties  were  examined  including  the  conduction  band,  valence  band, 
Fermi  level,  and  band  gap  energy.  The  difference  between  direct  and  indirect 
semiconductors  was  also  explained. 

Next,  the  PN-junction  was  introduced  opening  the  door  to  further  discussion  of 
semiconductor  devices  in  the  context  of  solar  cells,  and  Poisson’s  equation  was  used  to 
solve  for  the  electric  field,  potential,  and  depletion  region  width  in  the  junction.  The 
charge  distribution  within  the  device  was  analyzed  to  prepare  for  nonequilibrium  effects 
caused  by  solar  irradiance. 

The  important  concepts  of  solar  radiation  in  the  AM  1 .5  spectrum  were  discussed 
to  lay  the  framework  for  the  solar  illumination  analysis.  Following,  the  fundamental 
concept  of  carrier  generation  caused  by  solar  irradiance  was  presented,  along  with  the 
opposite  process  of  carrier  recombination. 

Finally,  the  solar  cell  photocurrent  was  investigated  to  create  the  total  current- 
voltage  relationship.  The  concepts  of  short  circuit  current,  open  circuit  voltage,  and  fill 
factor  were  defined  as  applied  to  the  characteristic  solar  cell  IV-curve.  To  wrap  up  the 
theoretical  analysis,  power  output  and  efficiency  were  described  and  defined  as  the  two 
most  important  measures  of  merit  for  practical  solar  power  systems. 

2,14  Literature  Review 

Since  the  goal  of  this  research  is  to  improve  the  modeling  of  solar  concentrator 
systems,  it  is  necessary  to  investigate  the  models  currently  used  in  industry  to  design 
these  systems. 

Following  the  familiar  equation  used  in  the  initial  analysis,  most  models  serve  to 
derive  both  components  of  the  current  using  the  assumption  of  constant  generation 
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throughout  the  device.  However,  this  is  an  oversimplification  of  the  total  current  since 
the  diode  current  and  photocurrent  are  inherently  dependent  upon  each  other  [14]. 
Furthennore,  unifonn  generation  throughout  the  device  is  not  a  practical  assumption, 
especially  given  the  wavelength  dependent  absorption  coefficients  for  most 
semiconductor  materials  compared  to  their  common  manufactured  thicknesses  for  use  as 
solar  cells. 

Simplified  assumptions  regarding  excess  carrier  generation  ultimately  lead  to  a 
simplified  expression  for  diode  current  throughout  the  device.  A  common  method  used  to 
determine  the  diode  current  in  a  PN-junction  solar  cell  is  to  solve  for  the  minority  carrier 
distributions  using  the  Shockley  method  for  a  short -base  //-side  and  long  base  //-side 
junction  without  incident  light,  shown  graphically  in  Figure  20  [12]. 


n-side  p-side 

P  n  n  p 


Figure  20:  Short  Base  (n-side)  Long  Base  (p-side)  PN-junction  [12] 

Referring  to  Figure  20  Wn  is  the  surface  of  the  device  (facing  the  sun),  x„  is  the 
depletion  region  boundary  on  the  n-side,  p’n  is  the  excess  concentration  of  holes  on  the  n- 
side,  n  ’p  is  the  excess  concentration  of  electrons  on  the  p-side,  -xp  is  the  depletion  region 
boundary  on  the  p-side,  -Lp  is  the  diffusion  length  for  electrons  on  the  p-side,  and -Wp  is 
the  physical  boundary  of  the  device.  The  model  shown  in  Figure  20  proves  valid  for  a 
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simple  semiconductor  PN-junctions  under  forward  bias  but  is  an  oversimplification  when 
used  for  solving  the  diode  current  in  a  solar  cell.  The  condition  of  photon  absorption  in 
solar  cells  introduces  excess  carriers  that  are  generated  in  the  bulk  of  the  surface  (n-side) 
layer.  Therefore,  the  excess  carrier  concentration  will  be  greater  than  zero  near  the 
surface  of  the  device  and  at  a  maximum  in  the  bulk  of  the  surface  layer.  The  resulting 
expression  for  diode  current,  derived  in  the  next  chapter,  is  much  more  complete  when 
surface  recombination  and  bulk  photon  absorption  is  included. 

Other  solar  cell  models  attempt  to  approximate  excess  carrier  concentration  in  a 
PN-junction  solar  cell  as  a  negative  exponential  throughout  the  device  as  shown  in  Figure 
21  [19]. 

n-side  p-side 


Figure  21:  Exponential  Model  for  Excess  Carrier  Concentration 

The  exponential  model  for  excess  carrier  concnetration  is  also  an 
oversimplification.  As  shown  previously  in  Figure  16,  photon  generation  can  be  modeled 
by  an  exponential,  but  the  phenonomenon  of  surface  recombination  (discussed  in  Section 
2.8)  serves  to  reduce  the  excess  carrier  concentration  at  the  surface  of  the  device  to  nearly 
zero.  Incorporating  the  effect  of  surface  recombination  dramatically  changes  the 
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resulting  equation  for  current  through  the  device,  as  will  be  shown  in  the  following 
analysis. 

Another  issue  with  modeling  excess  carrier  concentration  through  the  device  as  an 
exponential  arises  when  the  diffusion  length  approaches  the  photon  penetration  depth. 
This  anomaly  can  be  seen  in  the  equation  for  the  photogenerated  current  using  this 
method,  where  Ja  is  the  photocurrent,  q  is  the  elemental  charge,  A  is  the  surface  area  of 
the  cell,  a  is  the  inverse  of  penetration  depth,  Ga  is  the  generation  rate  at  the  surface,  and 
Ln  is  the  diffusion  length  for  holes  on  the  n-side  (absorption  layer)  [19]. 


2  1 

A  mathematical  singularity  arises  when  Ln  —  —  since  the  denominator  of 

equation  (2-15)  becomes  zero,  and  the  photogenerated  current  becomes  infinite.  The 
exponential  model  for  excess  carrier  concentration  was  initially  used  as  the  foundation 
for  this  research  until  this  singularity  was  discovered.  Other  authors  have  alluded  to  the 
relationship  between  diffusion  length  and  photon  penetration  depth,  but  it  was  confirmed 
that  when  mulitple  suns  are  used  to  illuminate  the  cell,  the  rise  in  cell  temperature  causes 
these  two  quantities  to  approach  equivalence  [22].  An  important  consequence  of  this 
discovery  is  the  realization  that  the  electrical  model  plays  an  important  role  in  detenning 
steady  state  temperature,  since  the  wrong  model  can  lead  to  drastically  extreme  values  for 
photogenerated  current  under  certain  conditions  leading  to  an  overestimation  of  the  heat 
generated  in  the  device.  Therefore,  further  investigation  into  the  excess  carrier 
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concentration  disrtibution  throughout  the  device  was  necessary,  and  a  new  model  was 
develped  to  serve  as  the  foundation  for  this  research. 

Several  simple  solar  cell  models  are  presented  in  the  academic  literature.  While 
some  of  these  models  oversimplify  significant  phenomena,  there  are  benefits  to  using 
such  simple  models  in  large-scale  applications.  In  situations  where  the  power  generated 
by  an  array  of  non-concentrated  terrestrial  solar  cells  exceeds  the  kilowatt  range,  absolute 
efficiency  becomes  less  important  than  raw  power  generated.  Accordingly,  when 
analyzing  these  systems,  elemental  dependence  upon  carrier  concentrations  under 
different  illumination  conditions  such  as  clouds  or  fog  proves  rather  insignificant  to  the 
broad  scale  application. 

Furthennore,  the  solar  cell  models  generated  using  the  solid  state  physics 
approach  cannot  accurately  predict  the  operation  of  non-crystalline  solar  cells. 

Amorphous  and  organic  cells  that  are  populating  the  research  as  third  generation  cells  are 
not  easily  predicted  using  detailed  numerical  analysis.  Analysis  of  these  cells  is  generally 
done  empirically  as  current-voltage  curves  can  be  drawn  for  any  photovoltaic  cell.  Also, 
the  physical  transition  from  low-level  injection  to  high-level  injection  as  a  function  of 
concentration  ratio  is  not  well  defined,  but  it  has  been  demonstrated  that  very  high  levels 
of  concentration  move  the  solar  cell  device  to  high-injection,  where  the  presented  model 
proves  invalid  [17]  [18]. 

The  most  significant  difference  between  the  solar  cell  models  currently  in  the 
literature  and  the  model  presented  in  this  research  is  that  the  model  presented  here 
quantifies  the  change  in  voltage  and  current  of  the  cell  as  a  function  of  the  dynamic 
thennal  effects  of  the  cell  under  illumination.  Under  concentration  the  cell  will  heat  up 
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and  the  device  parameters  will  change,  but  there  are  few  models  that  adequately 
investigate  this  phenomenon.  No  other  solar  cell  model  currently  in  the  literature  has 
solved  the  transport  and  continuity  equations,  simultaneously  with  the  heat  equations,  to 
investigate  the  effect  of  temperature  on  the  important  device  parameters  under  low-level 
injection.  Consequently,  this  research  is  beneficial  to  the  Air  Force  engineering 
community  especially  given  that  the  Air  Force  has  begun  installation  of  several  CSP 
systems  across  the  country. 

This  analysis  serves  to  forgo  the  simplification  used  often  in  photovoltaic  current 
solutions  by  solving  the  transport  and  continuity  equations  using  a  gamma  distribution 
approximation  for  the  excess  minority  carrier  concentration  throughout  the  device.  As 
will  be  seen,  the  total  current  derived  using  this  approach  is  not  nearly  as  aesthetic  or 
straightforward  to  derive,  but  as  with  most  complex  mathematical  models,  more 
accurately  predicts  the  physical  operation  of  the  device.  In  addition,  a  more  accurate 
current  model  lends  to  a  better  foundation  for  the  final  dynamic  thermal  model  of  the 
CSP  system,  as  will  be  shown  in  the  following  analysis. 


41 


TIT:  Gamma  Distribution  Excess  Minority  Carrier  Concentration  Model 

Accounting  for  exponentially  decaying  generation  using  a  gamma  distribution 
approximation  for  excess  minority  carrier  concentration  proves  complex.  However,  as 
will  be  shown,  the  resulting  differential  equations  can  be  analytically  solved  to  yield  a 
more  accurate  representation  of  the  minority  carrier  concentrations  and  resulting  current 
through  the  device. 

3. 1  Ideal  Diode  Analysis 

The  ideal  diode  equation  can  be  used  to  solve  for  the  total  current  through  the  PN- 
junction  solar  cell.  Solar  illumination  introduces  the  nonequilibrium  processes  of  carrier 
generation  and  recombination.  While  no  external  bias  is  applied  to  the  junction, 
generation  of  carriers  within  the  electric  field  of  the  depletion  region  effectively  causes 
the  PN-junction  to  be  forward  biased.  To  model  the  electrical  characteristics  of  a 
photovoltaic  cell,  a  solution  to  the  time  dependent  diffusion  equation  for  a  single  PN- 
junction  is  investigated.  The  device  is  analyzed  using  the  ideal  diode  equation.  The  goal 
of  this  analysis  is  to  characterize  the  solar  cell  by  investigating  the  photogenerated  diode 
currents  within  the  PN-junction  under  solar  illumination.  Since  ideal  diode  currents  are 
driven  by  minority  carriers,  the  excess  minority  carrier  concentration  on  each  side  of  the 
PN-junction  must  be  determined. 

3.2  Solar  Cell  PN-junction 

Generally,  a  simple  semiconductor  solar  cell  can  be  modeled  as  a  PN-junction 
with  a  relatively  large  depletion  region  and  an  asymmetric  base  length  for  the  n  and  p- 
side.  An  n-type  semiconductor  material  is  commonly  used  for  the  thinner  absorption 
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layer,  while  p-type  material  usually  comprises  the  thicker  base  region  because  the 
diffusion  length  for  electrons  is  much  longer  than  the  diffusion  length  for  holes.  Thus, 
creating  a  thicker  p-type  region  allows  for  maximum  diffusion  of  minority  carriers  [21]. 

Following  the  analysis  demonstrated  by  Muller  and  Kamins,  the  n-side  of  the 
junction  (facing  the  sun)  is  treated  as  a  short  base  case  since  the  diffusion  length  is  over 
the  entire  depth  of  the  junction  [14].  Accordingly,  the  p-side  is  treated  as  a  long  base 
case  since  the  diffusion  length  is  much  less  than  the  width  of  the  cell  [11].  Figure  22 
shows  a  diagram  of  the  cell  under  solar  illumination. 
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Figure  22:  PN-junction  Solar  Cell  Under  Illumination 


It  is  important  to  note  that,  given  the  condition  of  low  level  injection,  the  diode 
current  in  a  PN-junction  is  driven  by  the  excess  minority  carrier  concentration  (holes  on 
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the  n-side  and  electrons  on  the  p-side).  Thus  in  order  to  solve  for  the  current  across  the 
PN-junction  we  need  to  know  the  minority  carrier  concentration  on  each  side. 


3.3  Minority  Carrier  Concentrations 

The  excess  minority  carrier  concentration  on  the  n-side  of  an  ideal  short  base 
diode  under  forward  bias  can  be  approximated  by  a  linear  function,  provided  there  is  no 
carrier  generation  present.  However,  the  phenomenon  of  carrier  photo-generation 
coupled  with  surface  recombination  lends  the  excess  minority  carrier  concentration 
distribution  to  approximate  a  gamma  function  (with  shape  parameter  equal  to  two). 

The  gamma  distribution  PDF  is  shown  below  [25]: 


y(x)  = 


1 

bar(a ) 


_X 

xa~1e  b 


(3-D 


In  the  gamma  PDF  equation,  a  is  the  shape  parameter,  b  is  the  scale  parameter, 
and  T  is  the  gamma  function  term.  The  interesting  conclusion-  that  the  gamma 
distribution  can  be  used  to  model  excess  minority  carrier  concentration  in  a  solar  cell-  is  a 
unique  and  novel  contribution  of  this  research.  Several  other  excess  minority  carrier 
concentration  distributions  were  investigated,  but  no  other  distribution  was  able  to 
account  for  both  the  high  surface  recombination  rate  and  penetration  of  photons  into  the 
bulk  of  the  device,  and  it  has  been  demonstrated  that  surface  recombination  is  a 
particularly  important  phenomenon  in  PN-junction  silicon  solar  cells  [25].  The  gamma 
distribution  makes  sense  physically  because  the  excess  minority  carrier  concentration  at 
the  surface  of  the  device  (where  recombination  is  high)  is  negligible,  while  the 
concentration  peaks  at  a  depth  into  the  device  where  generation  is  high  and 
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recombination  is  low,  before  exponentially  decaying  throughout  the  rest  of  the  device.  A 
plot  of  the  gamma  distribution  reveals  the  physical  approximation  of  the  excess  minority 
carrier  concentration  on  the  n-side  of  a  solar  cell. 

On  the  /7-side  of  the  junction  the  excess  minority  carrier  concentration  can  be 
approximated  by  an  exponential  function  [11].  A  plot  of  the  excess  minority  carriers  for 
a  forward-biased  solar  cell  with  photo -generation  present  is  shown  in  Figure  23,  where 
p  and  n  ),  are  the  excess  holes  and  electrons  respectively,  x„  and  xp  are  the  depletion 
region  edges,  and  Wn  and  Wp  are  the  boundaries  of  the  device. 

n-side  p-side 


Figure  23:  Minority  Carrier  Concentration  Across  the  PN-junction 


The  plot  of  the  minority  carrier  concentrations  shown  in  Figure  23  assumes  bulk 
recombination  and  incorporates  exponential  EHP  generation,  with  the  sun  facing  the  n- 
side.  The  magnitude  of  the  generation  Go  caused  by  illumination  exponentially  decreases 
through  the  one-dimensional  depth  of  the  material  (x-direction  into  the  PN-junction),  as 
modeled  by  G0e~ax  and  shown  in  Figure  24. 
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Figure  24:  Exponential  Generation  [14] 


Furthermore,  while  no  external  applied  voltage  is  present  on  the  n-side,  the  solar 
illumination  that  generates  carriers  causes  the  PN-junction  to  be  forward  biased.  To 
investigate  the  flow  of  charge  through  the  device  under  forward  bias  given  the  excess 
minority  carrier  concentrations  shown  above,  it  is  necessary  to  solve  the  fundamental 
charge  transport  equations  on  both  sides  of  the  PN-junction  solar  cell. 

3.4  Continuity  Equation 

The  differential  continuity  equation  for  the  minority  carrier  concentration  on  the 
n-side  of  the  solar  cell  is  shown  below  [18]: 
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In  the  above  equation  p  ’n  is  the  minority  carrier  concentration  on  the  n-side 
(holes),  Dp  is  the  diffusion  coefficient  for  holes,  pp  is  the  mobility  for  holes,  E(x)  is  the 
electric  field  in  the  region,  g  ’(x)  is  the  generation  tenn  for  holes  in  the  n-side,  xpo  is  the 
carrier  lifetime,  and  x  represents  the  one-dimensional  distance  into  the  region. 

An  analogous  equation  can  be  written  for  minority  carriers  on  the  p-side  n  ’p  [20]. 
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Since  we  are  solving  for  the  steady  state  condition  of  the  solar  cell,  the  time  based 
derivative  of  minority  carrier  concentration  on  the  n-side  is  set  equal  to  zero. 

Furthennore,  using  the  depletion  approximation,  the  electric  field  in  both  the  n  and  p-side 
(outside  the  depletion  region)  can  be  set  to  zero. 
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Dividing  this  equation  by  the  diffusion  coefficient  Dp  yields  a  familiar  fonn  for  a 
second  order  homogenous  differential  equation,  where  the  diffusion  length  Lp  is  given  as 


'pXpo- 
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(3-5.1) 


d2(p'n)  g'(x )  p'n 

dx2  +  Dp  Lp2 

d2(n'p )  g'  (x)  n'p  (3-5.2) 

dx 2  +  Dn  Ln2 

As  shown  above  in  Figure  24  the  generation  tenn  g’(x)  can  be  modeled  as  a 
decaying  exponential  term.  The  generation  rate  will  be  highest  at  the  surface  of  the  solar 
cell,  since  the  light  intensity  will  be  greatest  at  the  surface,  and  the  rate  depends  upon  the 
absorption  coefficient  a  (material  property)  and  the  initial  generation  rate  Go  (function  of 
how  strong  the  light  is  at  the  surface).  Accordingly,  the  generation  rate  can  be  modeled 
by  the  tenn  g'  (x)  =  G0e~ax  [9]. 

The  resulting  differential  equation  must  be  solved  for  p  ’n  to  detennine  the 
minority  carrier  concentration  on  the  n-side  assuming  high  surface  recombination  and 
exponentially  decaying  generation  [13],  It  is  important  to  realize  that  the  entire  solution 
for  the  minority  carrier  concentration  on  the  n-side  will  be  the  sum  of  the  particular  and 
homogenous  solutions. 

3.5  Particular  Solution  to  the  Transport  Equation 

To  detennine  a  particular  solution  for  p  ’n  in  the  second  order  differential  equation 
shown  above,  a  conjecture  must  be  made  to  the  form  of  the  solution.  It  follows,  from  the 
form  of  the  second  order  differential  equation  shown  in  (3-5.1),  that  a  gamma  PDF 

n  p2 

p'n  =  ^—xe~8x  is  a  reasonable  guess  to  the  form  of  the  particular  solution  for  carrier 
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concentration,  where  5  and  [3  are  “dummy”  variables  used  to  satisfy  the  mathematical 
operations.  This  assumption  can  be  verified  by  plugging  the  particular  solution  into  the 
second  order  differential  equation  shown  in  equation  (3-6).  The  assumption  is  valid  as 
long  as  the  particular  solution  yields  a  non-trivial  solution  for  the  minority  carrier 
concentration. 

The  first  step  in  determining  the  validity  of  the  gamma  distribution  assumption  is 
to  solve  equation  (3-6)  using  the  minority  carrier  concentration  p  ’n  on  the  n-side.  The 
easiest  way  to  solve  this  equation  is  to  first  solve  for  each  of  the  three  parts  of  equation 
(3-6). 
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The  exponential  terms  can  be  eliminated  by  dividing  into  each  tenn  of  equation 
(3-10).  When  the  three  parts  of  equation  (3-6)  are  plugged  in  using  the  particular  solution 
as  shown  in  equation  (3-10),  a  non-trivial  solution  for  (3  (independent  of  x)  is  found  when 
5  =  a. 
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(3-11) 


In  the  second  order  differential  equation  for  excess  minority  carrier  concentration 

c  B82xe  ^ 

shown  above,  carrier  recombination  is  modeled  by  the  term - — ,  where  c0  —  — — —  is 

Tpo  r(2) 

the  recombination  coefficient.  Furthermore,  a  gamma  distribution  function  was  used  to 
solve  the  differential  equation  since  it  was  shown  that  surface  recombination  is  high  in 
PN-junction  solar  cells.  To  fundamentally  understand  the  mechanism  of  recombination 
in  terms  of  the  carrier  transport  equation  it  is  necessary  to  solve  for  the  recombination 
coefficient  c0  in  terms  of  photon  penetration. 

As  will  be  shown,  the  recombination  coefficient  invokes  a  correction  factor  that 
allows  a  solution  to  the  differential  equation  using  the  gamma  function  to  be  found. 
Without  the  recombination  coefficient  figured  into  the  transport  equation,  the  model  fails 
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to  satisfy  the  continuity  equation,  and  the  gamma  function  “guess”  does  not  solve  the 
second  order  differential  equation.  Therefore,  this  recombination  coefficient  is 
imminently  important  in  the  rest  of  the  model. 

To  investigate  the  recombination  term,  the  second  order  differential  terms 
dependent  upon  x  must  be  separated. 


(382(x82)  (382x 
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With  5  =  a,  an  expression  for  the  recombination  coefficient  cG 

in  terms  of  photon 

penetration  and  diffusion  length  can  be  found. 
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The  recombination  coefficient  shown  above  allows  the  transport  equations  to  be 
solved  using  the  gamma  distribution. 

3.7  Homogenous  Solution  the  Transport  Equation 

With  a  particular  solution  to  the  original  second  order  differential  equation,  the 
homogenous  solution  must  be  detennined  to  derive  an  expression  for  the  entire  minority 
carrier  concentration  on  the  n-side.  The  fonn  for  homogenous  solution  for  this 
differential  equation  is  shown  in  equation  (3-19)  [13]. 
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(3-19) 


—  (*-■ *n)  (*~*n) 

p'n  -  Ae  lp  +  Be  lp 
To  determine  the  homogenous  solution  for  the  minority  carrier  concentration  on 
the  n-side  it  is  necessary  to  use  boundary  conditions  to  solve  for  the  coefficients  A  and  B 
[11]. 


The  first  boundary  condition  used  to  solve  for  the  homogenous  solution  to  this 
second  order  differential  equation  is  derived  from  the  ideal  diode  law,  which  states  that 
no  excess  minority  carrier  concentration  is  present  at  the  outside  edge  xn  +  xp  (surface)  of 
the  material.  Therefore,  the  minority  carrier  concentration  at  the  surface  is  given  by 
equation  (3-20). 
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With  a  solution  to  one  of  the  coefficients,  a  second  boundary  condition  must  now 
be  applied.  The  second  boundary  condition  relies  on  the  depletion  approximation  and 
states  that  no  recombination  is  present  at  the  edge  of  the  depletion  region  xn. 
Furthennore,  at  the  depletion  region  boundary  the  excess  minority  carrier  concentration 

G„r(  2) 

due  to  generation  will  be - -  per  the  particular  solution,  but  the  total  excess  minority 

2  Dp  cc 


carrier  concentration  will  be  dependent  upon  the  magnitude  of  the  forward  bias  Va  caused 
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by  the  generation  of  carriers  under  illumination.  The  result  of  the  second  boundary 
condition  yields  another  solution  for  the  coefficients. 


p'n(xn)  =  pno 


( ekr  -  l 
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GJ\ 2) 
+  2  Dpcc3 
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Using  the  coefficient  equation  derived  using  the  first  boundary  condition,  a 


solution  for  both  coefficients  can  be  found. 


A  =  -Be  lp 


(3-25) 
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3.8  Total  Excess  Minority  Carrier  Concentration  on  the  n-side 


Once  the  two  coefficients  have  been  found,  a  final  expression  for  the  excess 

minority  carrier  concentration  on  the  n-side  as  a  function  of  one  dimensional  distance  x 

into  the  material  and  intensity  of  light  (quantified  by  Va)  can  be  defined. 

-(*-*")  (x~xn)  G  r(2)  (3-29) 

p'n(x)  —  Ae  lp  +  Be  lp  -\ - — 

H  2 D„a3 
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Now  that  the  current  on  the  n-side  of  the  solar  cell  PN-junction  has  been 
characterized,  the  current  on  the  p-side  must  be  investigated.  As  stated  earlier,  the 
continuity  equation  says  that  the  total  current  throughout  the  device  must  be  constant 
under  steady  state  conditions,  therefore  it  is  most  convenient  to  solve  for  the  two  currents 
at  the  depletion  region  boundaries  xn  and  -xp  and  sum  these  components  to  get  the  total 
current  through  the  device. 

3,10  Constant  p-side  Generation  Assumption 
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To  analyze  the  p-side  of  the  solar  cell  PN  junction  another  important  assumption 
is  made:  the  generation  rate  is  constant  throughout  the  />sidc  of  the  junction.  This 
assumption  proves  practical  given  the  physical  construction  of  most  semiconductor  solar 
cells,  where  two  differently  doped  layers  are  brought  into  contact  one  on  top  of  the  other 
as  shown  in  Figure  22.  Accordingly,  given  the  one-dimensionality  of  the  device,  most  all 
of  the  carrier  generation  will  occur  on  the  top  layer  of  the  PN  junction,  and  for  the  sake  of 
this  analysis  the  carrier  generation  throughout  the  bottom  p-side  layer  is  assumed  to  be  y. 

3,11  Continuity  Equation  for  p-side 

The  homogenous  equation  for  the  p-side  excess  minority  carrier  concentration  is 
rewritten  below: 

n  d2(n'p)  i  g'(x)  n'p  (3-33) 

°  =  ~d^  +  ~D^~L 7 

The  assumption  that  generation  is  constant  lends  a  much  more  aesthetic  solution 
to  this  second  order  differential  equation.  Furthermore,  the  solution  for  the  excess 
minority  carrier  concentration  on  the  long  base  p-side  (assuming  bulk  recombination  and 
constant  generation)  has  been  well  documented  in  the  literature  and  an  adequate  analysis 
was  followed  by  Honsberg.  The  resulting  general  solution  to  this  differential  equation  is 
shown  below  in  equation  (3-34),  where  y  denotes  the  excess  minority  carrier 

G  L  2 

concentration  due  to  constant  generation  0  71  . 

Dn 

( X+XP )  ~(x+xp)  (3-34) 

n'p  =  Ae  Ln  +Be  Ln  +y 
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g  r f2^ 

It  is  important  to  note  that  ^  J,  is  the  carrier  concentration  due  to  generation  at 

the  depletion  region  edge  of  the  n-side  of  the  junction,  and  the  assumption  follows  that 
the  carrier  concentration  due  to  generation  is  also  constant  throughout  the  depletion 
region  and  the  p-side  of  the  junction.  However,  the  form  of  the  carrier  concentration 
tenn  due  to  generation  will  change  slightly.  Following  a  similar  boundary  condition 
analysis  done  for  the  n-side  of  the  PN-j unction,  a  solution  for  the  minority  carrier 
concentration  as  a  function  of  x  on  the  p-side  is  derived.  The  first  step  is  to  solve  for  the 
minority  carrier  concentration  at  the  depletion  region  boundary  —xp  to  solve  for  the 
homogenous  coefficients  A  and  B  using  boundary  conditions  similar  to  ones  used  in  the 
analysis  on  the  n-side  of  the  junction.  The  first  boundary  condition  says  that  the  excess 
minority  carrier  concentration  will  be  greatest  at  the  depletion  region  edge  -xp  once  an 
external  voltage  (due  to  illumination)  is  applied.  Again,  y  refers  to  the  constant 
generation  throughout  the  bulk  of  the  p-side. 


,  (  £Ya 

)  =  npo  (^e  kr  - 

~{xb)  (xb) 

(3-35) 
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np0  [ekT  -  1 

)  -  Ae  lp  +  Be  lp 

(■ 2xb ) 

(  qva  \ 

(3-37) 

A  -  -Be  lp 

+  npo  lekT  —  1 1  e  lp 

The  second  boundary  condition  says  the  minority  carrier  concentration  on  the  p- 

side  at  the  far  end  of  the  junction  past  the  diffusion  length  (deep  into  the  device)  will  be 

near  the  equilibrium  concentration  npo  along  with  the  constant  generation  term  y. 

iM  (Zbl  (3-38) 

n'p(Lp)  —  y  —  Ae  lp  +  Be  Lv  +  y 
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e  kT 

V 

-  1  j  =  -Be  Lv  +  Be  Lv 

B  =  B  =  0 

(3-42) 

Therefore  B  must  equal  zero  to  satisfy  the  homogenous  solution  to  the  second 
order  differential  equation  shown  above,  which  follows  from  the  long  base  case 
coefficient  analysis  demonstrated  by  Muller  in  Kamins.  Furthermore,  upon  inspection 

-(x+Xp) 

this  seems  appropriate  because  of  the  exponential  term  Be  lp  following  the 
coefficient  B  in  the  homogenous  equation.  At  first  glance,  this  tenn  appears  to  be  a 
negative  exponential  but  reviewing  the  reference  dimensions  used  in  Figure  22  reveals 
that  indeed  this  tenn  will  be  a  positive  exponential  once  the  proper  negative  x-position 
into  the  material  is  used  in  the  solution.  It  follows  that  there  should  not  be  any  positive 
exponential  terms  in  the  excess  minority  carrier  concentration  equation  at  the  far  end  of 
the  diffusion  length  into  the  p-side  material. 

3,12  Excess  Minority  Carrier  Concentration  on  p-side 

With  both  coefficients  considered,  the  resulting  solution  to  the  second  order 
differential  equation  can  be  found. 
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It  has  been  shown  that  the  excess  minority  carrier  concentration  on  the  p-side  is 
also  a  function  of  the  magnitude  of  forward  bias  Va  caused  by  solar  illumination.  With  an 
expression  of  the  excess  minority  carrier  concentration  on  the  p-side  assuming  constant 
generation,  an  expression  for  the  current  density  through  the  p-side  of  the  junction  under 
forward  bias  caused  by  illumination  can  be  written. 


Jn  R^n 


dn'p 

dx 


Jn  —  R^n 


3.13  Total  Diode  Current 


(3-45) 

(3-46) 


The  diode  current  can  be  effectively  described  by  the  sum  of  the  currents  on  both 
the  n  and  p-side  of  the  junction,  at  each  depletion  region  boundary.  A  complete 
expression  for  the  current  density  through  the  device  that  accounts  for  the  contribution  of 


both  the  />sidc  and  n-side  is  written  below. 
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The  total  current  Itot  (A)  can  be  found  by  multiplying  the  current  density  by  the 
effective  area  of  the  solar  cell.  Again,  it  is  important  to  note  that  these  equations  hold  as 
long  as  the  device  operates  under  the  condition  of  low  level  injection.  The  dependence  of 
the  minority  carrier  concentration,  and  thus  total  current,  on  Va  uncovers  an  ostensibly 
contradictory  relationship  (since  we  detennined  that  p  was  not  a  function  of  applied 
voltage),  but  upon  further  investigation  an  interesting  phenomenon  is  explained.  With  no 
external  applied  voltage  the  diode  current  is  just  a  function  of  this  generation.  If  no 
generation  were  present,  as  is  the  case  with  the  classical  PN-junction  under  thennal 
equilibrium,  then  no  net  diode  current  would  be  present  without  an  applied  voltage.  This 
is  the  case  of  the  solar  cell  in  the  dark. 
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IV:  Thermal  Analysis 


As  the  solar  cell  collects  light  from  the  sun,  especially  under  concentration,  the 
physical  temperature  of  the  cell  rises.  Explained  earlier  in  the  analysis,  the  process  of 
EHP  generation  creates  excess  energy  that  is  expended  through  the  cell  as  heat.  Solving 
for  the  current  in  the  cell  allows  for  further  investigation  into  the  dynamic  thennal 
properties  of  the  semiconductor  solar  cell  under  concentration.  The  current  produced 
through  the  PN-junction  solar  cell  produces  heat  generated  due  to  solar  illumination  and 
the  associated  electrical  output.  Fundamentally,  from  an  electronics  perspective,  a  PN- 
junction  solar  cell  behaves  as  a  resistor  to  current  flow  when  an  external  circuit  is 
connected.  As  the  photocurrent  through  the  junction  increases  (illumination  and  carrier 
generation  increases)  more  power,  and  thus  heat,  is  created  throughout  the  device.  The 
total  current  provides  the  framework  from  which  heat  through  the  solar  cell  can  be 
investigated. 

4.1  Temperature  Effects 

With  an  expression  for  total  current  through  the  PN-junction  solar  cell,  accounting 
for  exponentially  decaying  generation  throughout  the  device,  thermal  properties  of  the 
cell  under  illumination  can  be  investigated.  Before  jumping  into  the  analytical  modeling 
of  heat  generated  throughout  the  solar  cell,  it  is  beneficial  to  understand  which  properties 
of  the  PN-junction  solar  cell  depend  upon  temperature  in  order  to  predict  what  properties 
are  going  to  change  once  the  temperature  of  the  cell  increases. 

Temperature  is  proportional  to  thennal  energy  within  a  semiconductor  [8].  As 
temperature  increases,  the  kinetic  energy  of  atoms  within  the  crystal  lattice  of  a 
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semiconductor  is  increased.  As  the  temperature  of  a  semiconductor  increases,  the  band 
gap  energy  Eg  decreases  [14].  Accordingly,  several  other  semiconductor  material 
parameters  are  affected  by  changes  in  temperature  as  well.  As  the  band  gap  energy  of  a 
semiconductor  decreases,  the  amount  of  photons  available  for  conduction  increases  since 
lower  energy  photons  can  then  be  readily  absorbed  under  increased  temperature.  It 
follows  that  lower  energies  are  required  to  break  the  covalent  bonds  within  the  material 
as  the  temperature  increases. 

An  increase  in  cell  temperature  causes  the  conductivity  a  to  increase  while 
decreasing  the  carrier  mobility  p,  which  increases  the  flow  of  current  within  a  PN- 
junction  solar  cell.  However,  this  increase  in  charge  flow  is  dominated  by  a  larger 
decrease  in  voltage  across  the  junction,  effectively  decreasing  the  overall  efficiency  of 
the  cell  as  the  temperature  increases.  Investigating  the  fundamental  device  equations 
quantitatively  describes  exactly  how  the  PN-junction  changes  with  temperature.  The 
fundamental  investigation  of  the  benefit  caused  by  increased  photon  generation  due  to 
solar  concentration  opposed  by  the  loss  of  efficiency  due  to  rising  cell  temperature  is  the 
crux  of  this  analysis.  As  will  be  shown,  adequate  cooling  of  concentrated  solar  power 
systems  is  vital  to  their  successful  operation. 


4.2  Diode  Current-Temperature  Relationship 

The  simplified  diode  current  Id  (derived  assuming  constant  generation)  contains 
an  observable  temperature  dependence  [12]. 
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As  the  temperature  of  the  cell  increases,  the  magnitude  of  the  diode  current  appears  to 
decrease  exponentially.  Further  investigation  of  the  equation  for  diode  current  reveals 
that  the  reverse  saturation  current  Iq  is  also  dependent  upon  temperature.  The  equation 
for  the  reverse  saturation  current  in  a  semiconductor  is  given  by  (4-2),  where  q  is  the 
elemental  charge,  A  is  the  area  of  the  cell,  Dp  and  Dn  are  the  diffusivities  of  holes  and 
electrons  (a  function  of  doping),  Lp  and  Ln  are  the  diffusion  lengths  holes  and  electrons, 
Na  and  No  are  the  doping  levels,  and  nt  is  the  intrinsic  carrier  concentration  (material 
property)  [12]. 


I0  =  qAni2 


Dn  Dp 

LnNo  LpNA 
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The  most  interesting  parameter,  for  the  sake  of  thennal  analysis,  in  the  equation 

for  reverse  saturation  current  (4-2)  is  the  intrinsic  carrier  concentration  n,  and  it  is 

assumed  in  the  equation  that  n,  »  No  and  NA.  The  intrinsic  carrier  concentration 

dependence  on  temperature  is  given  in  equation  (4-3),  where  m*e/h  are  the  effective 

electron  and  hole  masses,  and  B  is  a  temperature  independent  material  property  [20]. 

2nkT ^ 
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The  equation  for  n,  is  quite  complex,  but  it  is  evident  that  the  intrinsic  carrier 
concentration  is  dependent  upon  both  the  band  gap  energy  of  the  semiconductor  and  the 
temperature.  The  expanded  equation  for  intrinsic  carrier  concentration  can  be  plugged 
into  the  saturation  current  equation  to  reveal  the  important  temperature  dependence  on 
both  the  saturation  and  diode  currents. 
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Therefore,  by  inspection,  it  can  be  seen  that  as  the  temperature  increases  the  diode 
current  increases  as  well,  validating  the  previous  assertion.  However,  the  increase  in 
diode  current  with  temperature  fails  to  portray  the  complete  interaction;  the  junction 
voltage  must  now  be  investigated. 


4.3  Voltage-Temperature  Relationship 

The  electric  field  across  the  junction  is  also  affected  by  temperature.  As  the 
temperature  increases,  the  electric  field  experiences  a  significant  reduction  in  magnitude. 
Furthennore,  as  the  temperature  increases  the  depletion  region  width  W  decreases  as  well. 
The  equation  for  open  circuit  voltage  shows  an  obvious  temperature  dependence,  where 
Iph  is  the  simplified  current  due  to  illumination  and  I0  is  the  reverse  saturation  current 
[14]. 


Vqc  =  /cTZn(^) 
*0 


(4-6) 


A  relationship  showing  the  quantitative  open  circuit  voltage  difference  for  two 
different  temperatures  assuming  constant  bandgap  energy  can  also  be  derived  [15]. 

7V  (4-7) 


V0c2  =  V0Cl  (jf)  +  Eg(  1-jr) 


It  is  not  immediately  apparent,  but  when  realistic  values  for  temperature,  open 
circuit  voltage,  and  band  gap  energy  are  input  into  the  equation  above  the  resulting  open 
circuit  voltage  decreases  as  temperature  rises.  Using  an  initial  open  circuit  voltage  of  0.6 
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V,  an  increase  in  temperature  of  40  °C  from  300  K  yields  a  12%  decrease  in  the  resulting 
open  circuit  voltage  for  Silicon.  The  magnitude  of  this  open  circuit  voltage  decrease  far 
overshadows  the  small  increase  in  diode  current  effectively  decreasing  the  electrical 
output  of  the  cell.  A  qualitative  assertion,  based  on  the  preceding  analysis  of  the  ideal 
PN-junction  solar  cell,  reveals  that  a  necessary  goal  of  the  solar  cell  designer  should  be  to 
reduce  the  cell  operating  temperature  to  increase  both  output  and  efficiency. 

4.4  Heat  Generated 

There  are  three  components  contributing  to  the  heat  generated  in  a  PN-junction 
solar  cell.  The  first  component  of  heat  is  due  to  the  photons  of  light  below  the  band  gap 
energy  that  do  not  have  enough  energy  to  excite  electrons  up  to  the  conduction  band. 
These  low  energy  photons,  while  not  contributing  to  conduction,  still  contribute  to 
thennal  losses  in  the  cell  due  to  the  atom  interactions  within  the  lattice,  causing  weak 
atomic  vibrations  that  generate  heat. 

It  is  important  to  note  that  heat  generated  (and  dissipated)  can  be  expressed  in 
tenns  of  power,  and  for  the  sake  of  the  thermal  analysis  “power”  will  refer  to  heat.  To 
derive  the  expression  for  power  generated  due  to  photon  absorption  PT  below  the  bandgap 
the  reflection  coefficient  R  must  be  mentioned,  a  material  property  that  describes 
quantitatively  how  many  photons  reflect  off  the  surface  of  a  semiconductor  material  (see 
Section  2.7).  The  power  generated  depends  on  the  reflection  coefficient  since  any 
photons  that  reflect  off  the  solar  cell  surface  will  not  contribute  significantly  to  heat 
generated. 
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While  most  commercial  solar  cells  use  some  form  of  anti-reflection  (AR)  coating, 
the  assumption  is  made  that  the  test  cell  is  free  of  any  AR  material,  simplifying  the 
modeling  of  reflection  in  the  simulation.  Furthennore,  all  of  the  effective  wavelengths 
below  the  bandgap  must  be  accounted  for  and  given  the  diagram  of  the  AM  1.5  solar 
spectrum,  it  is  clear  that  there  will  be  a  large  range  of  wavelengths  emitted  by  the  sun 
below  the  band  gap  of  most  semiconductors.  To  include  all  effective  wavelengths,  the 
wavelength  dependent  generation  must  be  integrated  from  the  lowest  wavelength  up  to 
the  band  gap  energy.  The  resulting  expression  for  the  power  generated  by  photons  below 
the  band  gap  is  shown  in  equation  (4-8),  where  A  is  the  effective  area  of  the  cell,  G(A)  is 
the  wavelength  dependent  generation,  and  R  is  the  reflection  coefficient. 

rAEg  (4-8) 

Pt=A(1-R)\  G(A)dA 
Jo 

The  second  component  of  heat  in  the  PN-junction  solar  cell  is  generated  by  the 
flow  of  current  through  the  cell  given  by  Pj.  Incident  photons  with  energy  at  the  bandgap 
excite  electrons  from  the  valence  band  up  to  the  conduction  band.  This  conduction 
produces  the  current  derived  in  the  preceding  analysis.  The  power  generated  due  to 
illumination  is  assumed  to  be  constant  across  the  area  of  the  cell.  This  seems  a  valid 
assumption  given  that  the  physical  area  of  the  concentrated  solar  cell  is  small  so  that  the 
focusing  optics  can  be  configured  to  uniformly  illuminate  the  surface  [21].  Non-unifonn 
illumination  in  concentrated  solar  cells  was  investigated  by  Luque  et  ah,  and  it  was 
concluded  that  these  effects  are  incredibly  complex  to  model  [22]. 
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Given  unifonn  illumination  and  EHP  generation,  the  component  of  heat  due  to 
illumination  is  given  by  the  power  dissipated  across  the  junction,  where  Itot  is  the  total 
solar  cell  current,  and  Rceii  is  the  equivalent  resistance  of  the  solar  cell. 

Pi  =  ItotRceii  (4-9) 


The  final  component  of  heat  generated  in  the  PN-junction  solar  cell  is  caused  by 
incident  photons  with  energy  greater  than  the  bandgap  energy  of  the  semiconductor.  As 
explained  earlier,  these  photons  excite  electrons  past  the  conduction  band,  forcing 
electrons  to  relax  back  to  the  edge  of  the  conduction  band.  The  excess  amount  of  energy 
above  the  band  gap  that  a  photon  has  is  expended  through  the  cell  as  power  given  by  Pp. 

The  derivation  of  the  expression  for  power  generated  by  photons  above  the  band 
gap  is  very  similar  to  the  expression  for  power  generated  due  to  photons  below  the  band 
gap.  However,  the  quantum  efficiency  P  must  now  be  included  since  not  all  of  the 
photons  above  the  band  gap  will  be  absorbed  by  the  semiconductor.  Furthermore,  to 
account  for  all  of  the  effective  wavelengths  of  solar  radiation,  the  wavelength  dependent 
generation  term  must  be  integrated  over  all  wavelengths  above  the  associate  band  gap 
energy  wavelength.  The  resulting  term  for  the  power  generated  due  to  photons  above  the 
band  gap  is  shown  in  equation  (4-10),  where  P  is  the  quantum  efficiency,  hv  is  the  energy 
of  light,  Eg  is  the  band  gap  energy  and  G(X)  is  again  the  wavelength  dependent 
generation. 


Pp  —  A/3  I  (hv  -  Eg)G(X)dX 


(4-10) 


Eg 


The  three  terms  contributing  to  the  power  generated  in  a  PN-junction  solar  cell 
derived  above  provides  an  adequate  description  for  the  entire  heat  generated  in  the  cell. 
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4.5  Heat  Dissipated 

The  total  temperature  in  the  PN-junction  solar  cell  is  not  only  a  function  of  power 
generated  in  the  cell  but  also  a  function  of  the  power  dissipated  throughout  the  cell. 

There  are  two  primary  mechanisms  that  serve  to  dissipate  power  in  a  concentrated  solar 
cell.  Black  body  emission  Pbb  is  the  first  mechanism  that  removes  heat  from  a  cell,  since 
any  object  that  absorbs  electromagnetic  radiation  also  dissipates  a  fraction  of  that 
radiation  as  heat.  In  reality,  a  true  black  body  radiator  cannot  exist,  but  the  law  still 
applies  for  solar  cells  since  a  significant  amount  of  radiation  incident  upon  the  surface  of 
the  cell  is  absorbed.  Since  this  phenomenon  is  an  emission  versus  absorption,  the 
effective  change  in  temperature  of  the  cell  is  negative. 

The  black  body  heat  phenomenon  was  first  described  by  Max  Planck,  and  an 
eloquent  solution  is  given  for  the  heat  generated  by  a  black  body  emitter  as  a  function  of 
temperature,  surface  area,  and  the  Stefan-Boltzmann  constant  a  [23]. 

Pbb=AoT 4  (4-11) 

The  second  mechanism  that  serves  to  remove  power  from  the  solar  cell  is  due  to 
convection  currents  created  by  active  cooling  systems.  These  active  cooling  systems  are 
standard  in  most  concentrated  solar  cells  and  are  found  in  a  variety  of  fonns  using 
various  coolant  liquids  such  as  air,  water,  and  other  fluids  [21].  For  the  sake  of  this 
analysis,  air  is  the  primary  coolant  and  thus  the  heat  convection  is  a  function  of  the  air 
temperature  and  flow  rate  across  the  surface  of  the  cell.  In  thermodynamic  analysis  of 
the  solar  cell,  the  physical  surface  of  the  cell  is  treated  as  a  simple  2-dimensional  plate. 
The  power  transferred  Pc  from  the  plate  to  the  laminar  air  stream  provided  by  the  active 
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cooling  system  is  given  by  the  expression  below,  where  hpiate  is  the  heat  transfer 
coefficient  of  the  semiconductor  material,  A  is  the  effective  area  of  the  cell,  Tceu  is  the 
temperature  of  the  cell,  and  Tm  is  the  temperature  of  the  laminar  cooling  stream  [23]. 

Pc  =  hpiate A(TCeU  ~  Poo)  (4-12) 

Another  heat  dissipation  method  used  in  most  solar  cells  is  the  physical  heat  sink. 
Using  a  heat  sink  to  dissipate  heat  from  the  cell  is  an  effective  measure  in  cells  with  large 
surface  areas.  However,  heat  sinks  must  also  dissipate  heat  through  some  form  of  a 
convective  cooling  mechanism.  It  is  also  important  to  note  that  heat  transfer  by  means  of 
conduction  was  neglected  from  the  total  heat  equation,  since  the  physical  width  of  the 
device  is  so  small. 

4.6  Total  Heat 

With  the  five  main  components  of  heat  in  the  solar  cell  accounted,  an  expression 
for  the  total  power  can  be  written. 

PT  =  PT  +  P,  +  Pp-  (Pc  +  PBB)  (4-13) 

When  the  total  power  approaches  zero  a  steady  state  temperature  point  is  reached 
where  the  power  generated  in  the  cell  is  exactly  balanced  by  the  power  dissipated. 
Initially,  the  power  generated  in  the  cell  will  be  greater  than  the  power  dissipated  but  as 
time  continues  and  all  other  factors  are  controlled  the  two  mechanisms  removing  heat 
from  the  cell  will  eventually  settle  the  solar  cell  temperature  to  a  new  equilibrium 
temperature.  It  is  important  to  realize  that  the  convective  cooling  system  serves  to  keeps 
the  solar  cell  from  heating  up  to  thermal  runaway,  a  disastrous  positive  feedback 
condition  for  semiconductors  in  which  the  cell  heats  up  to  failure  [25].  As  long  as  the 
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cooling  system  is  sufficient,  the  convective  cooling  apparatus  will  remove  a  significant 
portion  of  the  heat  generated  and  a  steady  state  temperature  will  be  reached. 

The  goal  of  the  dynamic  thennal  model  is  to  iteratively  solve  the  continuity 
equation  for  different  temperatures,  accounting  for  the  five  components  of  power 
discussed  above,  in  order  to  analytically  determine  exactly  when  this  steady  state 
temperature  point  is  realized. 

The  value  of  this  analysis  is  that  it  can  quantitatively  detennine  optimum 
concentration  factors  for  a  CSP  system  under  low  levels  of  concentration  because  the 
only  factor  that  will  fundamentally  change  the  dynamic  nature  of  the  model  (and 
effective  increase  in  temperature)  is  the  magnitude  of  solar  irradiation.  Since  it  was 
detennined  previously  in  the  analysis  that  the  photocurrent  generated  in  a  cell  linearly 
increases  with  solar  irradiation,  while  the  effective  increase  in  temperature  serves  to 
decrease  the  electrical  output,  an  optimum  balance  between  concentration  ratio  and 
temperature  must  be  achieved  assuming  that  the  active  cooling  mechanism  remains 
constant.  The  thennal  model  will  be  run  for  varying  levels  of  concentration  to  investigate 
the  associated  increases  in  cell  temperature.  Obviously,  lower  concentration  factors  will 
lead  to  lower  cell  temperatures,  and  the  thermal  model  will  reach  steady  state  more 
quickly  than  high  concentration  ratios.  Additionally,  at  lower  concentrations  the 
magnitude  of  the  current  will  be  lower  than  for  higher  concentrations.  The  goal  is  to 
detennine  where  the  two  quantities  are  optimized. 

Another  important  consideration  stems  from  a  fundamental  assumption  used  in 
this  analysis,  the  assumption  of  low-level  injection  operating  conditions.  Under  low  level 
injection,  changes  in  majority  carrier  concentrations  are  insignificant  relative  to  the 
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minority  carrier  concentrations.  Thus,  the  majority  carrier  concentrations  on  either  side 
of  the  PN-junction  remain  constant.  However,  the  process  of  concentration  serves  to 
drastically  change  the  carrier  concentrations,  as  shown  previously,  and  at  very  high  levels 
of  concentration  low  level  injection  no  longer  applies.  Once  the  device  changes  from  low 
to  high  level  injection  conditions,  most  all  of  the  preceding  analytical  derivations  for 
carrier  concentrations,  depletion  region  width,  and  electric  field  prove  invalid. 

Therefore,  this  thermal  model  only  deals  with  low  level  injection  and  accordingly  only 
low  concentration  ratios. 

4.7  Thermal  Model 

The  fundamental  idea  behind  the  dynamic  thermal  model  of  a  PN-junction  solar 
cell  under  concentration  is  relatively  straightforward.  The  continuity  equation  derived  a 
PN-junction  under  illumination  is  solved  under  steady  state  conditions  for  one 
concentration  and  a  range  of  different  temperatures,  and  the  total  power  is  calculated  at 
each  temperature.  The  temperature  at  which  the  power  generated  equals  the  power 
dissipated  (thus  total  power  equals  zero)  is  the  steady  state  temperature.  Furthermore, 
IV-curves  are  plotted  for  each  temperature  to  easily  investigate  how  the  fill  factor  and 
open  circuit  voltage  change  with  temperature.  The  MATLAB  code  generated  for  the 
thennal  model  is  included  in  Appendix  A. 
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V:  Test  Cell  Characterization 


In  order  to  verify  the  accuracy  of  the  theoretical  thermal  model  presented  in  this 
work,  it  is  necessary  to  conduct  laboratory  tests  on  a  solar  cell  under  maximally 
controlled  conditions.  Furthennore,  to  ensure  that  the  conducted  tests  are  as  thorough  as 
possible,  it  is  necessary  to  completely  characterize  the  material  properties  of  the  test  cell 
to  input  into  the  computer  model,  so  that  the  two  results  (physical  and  theoretical)  can  be 
simultaneously  compared  and  analyzed. 

The  solar  cell  used  throughout  testing  is  a  simple  commercially  available 
polysilicon  cell.  A  chronological  and  reproducible  series  of  steps  was  performed  to 
characterize  the  material  properties  of  the  test  cell.  First,  the  junction  capacitance  of  the 
cell  was  measured  to  investigate  the  impurity  doping  concentrations  of  the  PN-junction. 
Next,  the  effective  cell  thickness  was  measured  using  a  scanning  electron  microscope 
(SEM).  Following  the  SEM  measurement,  the  reverse  saturation  current  was  calculated 
from  Na  and  N(|  to  investigate  the  n-side  and  p-side  thicknesses  individually.  The  final 
characterization  step  was  perfonned  using  the  doping  concentrations,  thicknesses,  and 
resistivities  to  detennine  the  internal  resistance  of  the  cell.  It  is  important  to  realize  that 
all  characterization  tests  of  the  cell  were  perfonned  with  the  cell  in  the  dark  to  maximize 
the  control  of  the  tests  and  to  evaluate  the  test  cell  as  an  ideal  diode. 

5,1  Phase  Shift  Method  to  Determine  Junction  Capacitance 

A  very  clever  method  for  detennining  doping  levels  indirectly  through  variations 
injunction  capacitance  was  demonstrated  by  Timken,  Milsted,  et  al.  A  diagram  of  the 
test  set-up  is  shown  in  Figure  25. 
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Figure  25:  CV  Measurement  Test  Diagram 

The  capacitance-voltage  (CV)  measurement  technique  is  an  established  technique 
to  accurately  determine  impurity  concentrations  in  PN-junction  solar  cells  [27],  and  this 
test  is  primarily  conducted  in  the  time  domain  by  measuring  the  transient  decay  of  a 
voltage  step  through  the  cell.  However,  this  method  proved  difficult  in  the  lab  because  of 
the  relatively  large  capacitance  of  the  voltage  source  used  to  generate  the  applied  voltage. 
Therefore,  this  CV-test  was  performed  in  the  frequency  domain  to  measure  the  phase 
shift  of  the  output  sinusoidal  waveform  from  the  solar  cell  (in  the  dark)  as  compared  to 
the  input. 

A  PN-junction  under  thermal  equilibrium  can  be  modeled  as  a  capacitor,  where 
the  parallel  plates  are  the  quasi-neutral  n  and  p-sides  of  the  device  separated  physically 
by  the  width  of  the  depletion  region.  The  theory  behind  the  test  is  grounded  in  the  fact 
that  the  capacitance  of  a  PN-junction  diode  changes  with  applied  voltage,  since  the  width 
of  the  depletion  region  also  changes  with  applied  voltage  as  discussed  in  Chapter  III. 
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qesNaNd 


(5-1) 


C'  = 


In  the  equation  above  VR  is  the  applied  voltage  across  the  cell.  The  capacitance 
of  the  cell  under  zero  applied  voltage  can  be  directly  measured  with  a  multimeter.  But, 
the  capacitance  of  the  cell  under  applied  bias  must  be  measured  indirectly  with  the 
frequency  domain  CV-test.  Using  the  set  up  shown  in  Figure  25,  the  phase  shift  of  the 
output  waveform  compared  to  the  input  was  measured  for  the  solar  cell  with  an  applied 
voltage  of  approximately  -6.6  V.  A  diagram  of  the  output  waveform  from  the 
oscilloscope  is  shown  in  Figure  26. 
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Figure  26:  Oscilloscope  Output  Waveform 

The  phase  difference  Arp  of  2.2  us  was  converted  to  a  capacitance  value  by 
equation  (5-2)  [24]: 
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(5-2) 


c  _  sin(A (p)  Vres 
2uf  Rm^cell 

Referring  to  equation  (5-2)  above,  vres  is  the  magnitude  of  the  voltage  across  the 
series  resistor,  Rm  is  the  series  resistance  value,  and  vcell  is  the  magnitude  voltage  across 
the  solar  cell.  With  a  forward  bias  of  -6.6  V  and  a  phase  difference  of  2.2  ps,  the 
capacitance  of  the  solar  cell  was  measured  to  be  174.3  nF,  a  reasonable  value  for  diode 
capacitance.  To  verify  the  accuracy  and  validity  of  this  test,  the  same  measurements 
were  taken  with  zero  applied  voltage,  and  the  capacitance  measured  from  the  phase  shift 
method  was  identical  to  the  multimeter  value  of  500.4  nF,  adding  faith  and  reliability  to 
this  method. 

With  two  capacitance  values  at  two  different  applied  voltages,  the  doping 
concentrations  could  be  investigated. 

5.2  Donor  and  Acceptor  Doping  Concentrations 

The  equation  for  capacitance  of  a  PN-junction  reveals  an  obvious  dependence 
upon  doping  [12]. 


qesNaNd 

A 

2 In  (“pr')  +  VR)  (Na  -1-  Nd) 

(5-3) 


This  relationship  can  also  be  expressed  in  tenns  of  products  P  and  sums  S  of 
doping  concentrations. 
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(5-4) 
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The  latter  form  lends  itself  well  to  numerical  analysis,  given  the  preceding  CV- 
test  measurements.  Since  two  capacitance  values  were  measured  at  two  different  values 
of  Vr,  a  system  of  two  equations  with  two  unknowns  can  be  solved  for  P  and  S. 

(5-5) 

(5-6) 

The  system  of  non-linear  equations  was  solved  directly  in  MATLAB  for  P  and  S 
yielding  a  value  of  approximately  2  x  1032  for  P  and  3.2  x  10 16  for  S'.  Given  that  the 
product  of  the  doping  values  was  so  much  larger  than  the  sum,  it  was  inferred  that  the 
impurity  doping  of  the  test  cell  was  approximately  equal  on  each  side  of  the  PN -junction, 
with  an  approximate  magnitude  of  1.7  x  10  cm'  . 
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applied  voltage  across  the  cell,  along  with  the  total  current  through  the  cell  by  the 
following  equation: 


It  kT 
14  =  In  (/+!)■ 


7C 


R 


(5-8) 


In  the  equation  above,  It  is  the  total  current  through  the  cell  and  Va  is  the  external 
voltage  applied  across  the  cell.  Ideally,  Is  remains  constant  with  changes  in  applied 
voltage  and  current;  therefore,  if  voltage  and  current  is  known  the  reverse  saturation 
current  can  also  be  found.  A  simple  circuit  was  devised  to  evaluate  this  test  condition,  in 
which  the  solar  cell  was  placed  in  series  with  a  5  V  DC  source  and  variable  resistor  Rl. 
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Figure  27:  Reverse  Saturation  Current  Test  Diagram 

The  total  current  through  and  voltage  across  the  cell  is  easily  measured  with  a 
multimeter  since  the  cell  is  in  series  with  everything  else  in  the  circuit.  However,  there 
are  two  unknowns  in  the  circuit  shown:  internal  resistance  rj  and  the  voltage  drop  only 
across  the  ideal  diode.  The  voltage  v;  measured  with  the  multimeter  includes  the  voltage 
drop  across  the  internal  resistance,  but  in  order  to  fully  investigate  the  physical  properties 
of  the  cell,  especially  the  reverse  saturation  current,  the  voltage  va  must  be  determined 
since  there  is  a  observable  relationship  between  va  and  Is. 
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Vi  =  val  +  hrd 


(5-9) 


u  kT 
=  In  (7"  +  1)  — 


yal 


(5-10) 


Since  there  are  two  unknown  values  Is  and  rj,  the  current  and  voltage  in  the  circuit 
must  be  measured  twice,  for  two  different  applied  voltages,  so  that  two  equations  with 
two  unknowns  can  be  solved. 


u  kT 
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(5-11) 

(5-12) 


The  reverse  saturation  current  tenn  must  be  iteratively  solved  using  the  Monte 
Carlo  method  for  two  different  applied  voltages  and  corresponding  currents  to  minimize 
the  measurement  error.  Using  the  Monte  Carlo  method  and  the  test  circuit  described,  the 
reverse  saturation  current  density  was  detennined  to  be  approximately  6.269  pA/cm“. 

The  value  for  resistance  rd  from  this  method  proves  incomplete,  however.  What  the 
quantity  rd  physically  represents  is  the  internal  resistance  of  the  solar  cell,  but  the  value  rd 
measured  in  the  preceding  analysis  includes  the  resistance  of  the  wires  and  contacts,  a 
significant  contribution  to  error  when  typical  diode  resistance  values  are  considered. 
Therefore,  the  precise  value  of  diode  resistance  must  be  detennined  another  way. 


5.4  N/P-side  Thickness  Determination 

In  order  to  evaluate  the  test  solar  cell  internal  resistance,  the  n-side  and  p-side 
thicknesses  of  the  PN-junction  must  be  further  investigated.  The  total  thickness  of  the 
cell  was  measured  using  the  scanning  electron  microscope  (SEM). 
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Figure  28:  SEM  Image  of  the  Test  Cell 

Since  the  reverse  saturation  current  density  of  the  cell  was  solved  in  the  previous 
characterization,  the  width  of  the  n-side  Wb  can  be  explicitly  solved. 
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Using  known  values  characterized  in  the  previous  analysis,  the  width  of  the  n- 
side  of  the  test  solar  cell  was  determined  to  be  40.7  pm.  Given  the  n-side  thickness, 
depletion  region  thickness  (xb  +  |xp|),  and  total  cell  thickness,  the  p-side  thickness  can 
easily  be  solved.  The  p-side  thickness  of  the  test  solar  cell  was  determined  to  be 
435.6  pm. 


5.5  Solar  Cell  Internal  Resistance 

The  internal  resistance  rd  of  a  PN-junction  diode  is  described  by  the  following 
equation  [12]: 


rd  =  j[pnthn  +  Ppthp\ 


(5-14) 
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Equation  (5-14)  is  an  applied  fonn  of  the  familiar  equation  for  resistance  The 

n-side  and  p-side  thicknesses  are  described  by  thn  and  thp  respectively,  while  p„  and  pp 
describe  the  resistivities  (as  a  function  of  doping)  of  both  regions  of  the  cell.  The  internal 
resistance  of  an  ideal  diode  is  usually  so  low  that  it  cannot  be  measured  with  a 
multimeter.  Therefore,  the  only  way  to  determine  diode  resistance  is  to  solve  the 
equation  above.  The  depletion  region  can  be  neglected  in  the  calculations  because  it  is  so 
small  relative  to  the  quasi-neutral  regions  of  the  solar  cell.  The  internal  resistance  of  the 
solar  cell  used  during  testing  was  detennined  to  be  approximately  2.8  mQ. 

5.6  Illumination  Test  Calibration 

A  tungsten-filament  halogen  lamp  was  used  for  the  duration  of  illumination 
testing.  An  important  issue  encountered  during  any  indoor  solar  cell  test  is  the  irradiance 
spectrum  mismatch  between  the  manufactured  light  source  and  the  sun.  A  comparison  of 
the  AM  1.5  spectrum  and  the  tungsten  lamp  spectrum  is  shown  in  Figure  29.  Using  the 
reference  cell  to  calibrate  the  light  source,  the  known  short  circuit  current  was  reached  for 
1  sun.  Under  this  lighting  condition,  a  short  circuit  current  of  253.9  mA  was  measured 
for  the  test  cell. 
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Figure  29:  Spectra  Mismatch  [26] 

No  light  source  available  exactly  matches  the  spectral  irradiance  pattern  of  the 
sun,  so  in  order  to  maximize  the  accuracy  of  the  dynamic  thermal  model  a  scaling  factor 
must  be  found  that  relates  the  light  source  spectrum  to  the  ideal  solar  irradiance  spectrum. 
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VI:  Illumination  Testing 


After  characterizing  the  physical  properties  of  the  test  solar  cell,  the  steady  state 
temperature  under  various  concentrations  was  measured. 

6.1  Test  Set-up 

The  cell  was  illuminated  using  150  Watt  tungsten  halogen  lamps.  A  small  electric 
fan  provided  a  laminar  flow  cooling  stream  across  the  surface  of  the  cell.  The 
temperature  of  the  cell  was  measured  using  a  thermal  couple  soldered  to  the  cell.  A 
second  multimeter  was  used  to  measure  the  short  circuit  current  to  detennine  the 
concentration  ratio  incident  on  the  test  cell,  given  that  short  circuit  current  is  linearly 
dependent  on  concentration  and  changes  relatively  minimally  with  increasing 
temperature.  A  diagram  of  the  test  set-up  is  shown  in  Figure  30. 


Figure  30:  Illumination  Test  Set-up 

However,  before  the  testing  could  begin  the  thennal  plate  convective  cooling 
properties  of  the  cell  were  measured  to  determine  the  value  for  hpiate,  to  include  in  the 
computer  model. 
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6.2  Cooling  Coefficient  Calibration 

The  convective  cooling  coefficient  tenn  hpiate  was  introduced  in  the  previous 
chapter,  and  is  a  fundamental  component  of  the  thermal  model.  Since  the  convective 
cooling  properties  depend  on  the  airflow  rate  over  the  solar  cell  surface,  hpiate  must  be 
calculated  for  the  laboratory  test  condition  fan  speed  in  order  to  maximize  control  in  the 
experiment.  The  convective  cooling  coefficient  given  in  equation  (6-1)  is  a  function  of 
the  specific  heat  of  the  solar  cell  material  hsu  (Silicon)  and  the  measured  exponential 
cooling  rate  of  the  cell  r/,.  An  extensive  derivation  of  equation  (6-1)  is  given  in  Appendix 
B. 

_hsu  (6-1) 

^ plate  ~ 

Th 

To  investigate  the  parameter  hpiate,  the  cell  was  heated  up  to  above  the  steady  state 
temperature  for  one  sun,  and  then  the  light  source  was  turned  off.  The  temperature  of  the 
cell  was  recorded  with  time  as  soon  as  the  light  source  was  turned  off,  essentially 
measuring  the  time  back  to  room  temperature  for  the  cell.  Plotting  the  temperature  of  the 
cell  as  a  function  of  time  the  exponential  form  of  convective  cooling  is  revealed  and  can 
be  modeled  by  equation  (6-2). 

-A  (6-2) 

T(t)  =  Tne 

In  equation  (6-2)  T(t)  is  the  temperature  as  a  function  of  time  in  the  convective 
cooling  experiment,  Tn  is  the  normalized  maximum  temperature,  t  is  time,  and  i|,  is  the 
cooling  rate.  The  temperature  was  plotted  as  a  function  of  time  for  the  maximum  fan 
speed  in  order  to  measure  Th. 
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Convecti\e  Cooling  Calibration  Test 


Figure  31:  Convective  Cooling  Calibration 

Referring  to  the  graph  shown  in  Figure  31,  the  r/,  is  6.623  from  the  equation  of  the 
exponential  regression.  The  convective  cooling  coefficient  hpiate  can  then  be  found  using 
equation  (6-1),  and  the  convective  cooling  part  of  the  model  is  complete.  As  an 
addendum,  r/,  for  five  different  fan  speeds  in  measured  in  Appendix  C.  Furthermore,  a 
trend  plot  of  hpiate  as  a  function  of  fan  speed  is  given  in  order  to  investigate  the 
importance  of  flow  rate  to  the  convective  cooling  properties  of  the  cell. 

A  MATLAB  plot  shows  the  extrapolated  data  points  demonstrating  higher  fan 
speeds  than  were  tested  in  the  lab  in  Figure  32. 
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Laminar  Air  Stream  Velocity  vs  Cooling  Parameter  th 


Figure  32:  MATLAB  Plot  of  Extrapolated  Cooling  Parameter  Values 

Obviously,  as  the  fan  speed  increases  the  cooling  coefficient  hpiate  increases  (eq. 
6-1,  6-2)  along  with  the  effective  cooling  of  the  cell.  However,  a  higher  fan  speed 
requires  more  power  in  order  to  run  the  fan,  resulting  in  an  optimization  problem 
concerning  both  cooling  and  power  requirements.  This  issue  presents  an  opportunity  for 
future  research  in  this  area.  After  the  convective  cooling  coefficient  of  the  cell  was 
characterized,  testing  under  various  concentration  ratios  commenced. 

6.3  Test  Procedures 

The  illumination  testing  procedures  were  relatively  straightforward.  There  were 
three  steps  used  throughout  testing  to  repeatedly  measure  steady  state  temperature  of  the 
cell  under  different  illumination  conditions: 

1.  Measure  the  short  circuit  current  (Isc)  at  one  sun 

2.  Use  the  measured  Isc  to  determine  different  concentration  ratios;  increase  the 
intensity  of  the  light  sources  accordingly  to  reach  desired  concentration  conditions 
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3.  Record  the  steady  state  temperature  for  different  concentration  ratios 


Following  the  three  steps  described,  the  steady  state  temperature  of  the  test  cell  was 
measured  as  a  function  of  the  concentration  ratio. 

6.4  Steady  State  Temperature  Measurements 

Figure  33  shows  the  theoretically  predicted  and  experimentally  measured  steady 
state  cell  temperature  as  a  function  of  concentration  ratio.  A  strong  agreement  is 
observed  between  the  two,  where  the  exponential  concentration  ratio  was  limited  to  four 
suns  due  to  limitations  of  the  cell  under  test. 


Experimental  and  Simulated  Results  Comparison 


Figure  33:  Illumination  Test 

The  maximum  steady  state  cell  temperature  measured  was  398  K  (257  °F),  at  a 
concentration  factor  of  four  suns.  Any  further  increase  in  light  intensity  resulted  in 
temperatures  that  compromised  the  test  set-up  and  damaged  the  test-cell. 
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VII:  Simulation  Results  and  Analysis 


7.1  Simulation  Predictions 

To  predict  the  steady  state  temperature  of  the  cell  as  a  function  of  concentration, 
the  thermal  model  was  run  for  different  values  of  concentration  and  the  temperature  was 
calculated.  Figure  34  shows  a  MATLAB  plot  of  steady  state  temperature  predictions  as  a 
function  of  concentration  up  to  35  suns. 


Cell  Temperature  as  a  Function  of  Concentration 


Figure  34:  Predicted  Cell  Temperature  as  a  Function  of  Concentration 

It  is  interesting  to  note  that  at  a  concentration  of  35  suns  and  temperature  of  720  K 
the  intrinsic  carrier  concentration  /?,  for  silicon  is  3.8  x  10  /cm" ,  approaching  the 
impurity  doping  levels  measured  from  the  test  cell:  Na  and  No  ~  1.7  x  1017.  From  a 
purely  theoretical  perspective,  the  simulated  cell  appears  to  violate  the  condition  of  low- 
level  injection  near  30  suns  of  concentration. 
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7.2  Effect  of  Temperature  on  the  Electrical  Properties 

The  thermal  model  developed  in  MATLAB  provides  a  tangible  demonstration  of 
the  severe  effects  an  increase  in  temperature  has  on  solar  cell  operation.  As  the  solar  cell 
heats  up  under  concentration,  the  open-circuit  voltage  suffers  while  the  short  circuit 
current  increases  linearly,  effectively  changing  the  shape  of  the  IV-curve  as  the 
concentration  is  increased.  Figure  35  shows  modeled  IV-curves  for  various  concentration 
ratios,  where  the  theoretically  predicted  steady  state  cell  temperature  is  incorporated.  For 
each  case,  an  impedance  matched  load  is  assumed,  maximizing  the  fill  factor  and  power 
delivered  by  the  cell. 


Simulated  IV-cunes  Under  Different  Concentrations 


Figure  35:  IV  Curves  for  5  Different  Concentration  Ratios 


First  order  models  for  concentrated  solar  cells,  not  accounting  for  changes  in 
temperature,  presume  a  simple  linear  relationship  between  concentration  ratio  and  power 
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output.  However,  as  shown  in  Figure  35,  changes  in  temperature  have  drastic  effects  on 
the  maximum  power  output  of  the  cell.  Ideally,  a  CSP  system  should  have  a  sufficient 
cooling  system  to  keep  the  open-circuit  voltage  high  even  under  high  concentration 
ratios.  However,  the  lab  set-up  demonstrates  that  the  cooling  system  used  was  largely 
ineffective  at  removing  enough  heat  to  keep  the  open  circuit  voltage  from  collapsing 
under  concentration.  The  fundamental  dependence  of  the  open  circuit  voltage  upon 
temperature  presents  an  important  future  research  opportunity  in  better  solar  cell  cooling 
methods. 

7.3  Optimum  Concentration  Factor  Analysis 

Using  the  multiparameter  approach  described,  the  optimum  concentration  factor 
of  a  CSP  system  by  incorporating  the  thennal  effects  of  increased  illumination  into  the 
electrical  model  of  solar  cell  operation.  The  predicted  optimum  concentration  for  the  cell 
used  in  the  test  set-up  is  plotted  in  Figure  36,  where  it  was  determined  that  2.9  suns  was 
the  optimum  concentration  factor  for  the  solar  cell  used  in  the  test-set  up. 
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IV-curve  for  Concentration  Factor  of  2.9  Suns 


Figure  36:  IY-curve  for  Optimum  Concentration 

Given  the  cooling  system  and  cell  material,  any  higher  concentration  ratios  would 
degrade  the  output  of  the  cell  due  to  the  collapse  of  the  open  circuit  voltage,  despite  the 
increase  in  light  intensity.  This  is  a  very  interesting  conclusion  because  most  commercial 
CSP  systems  use  concentration  ratios  on  the  order  of  hundreds  of  suns.  Therefore,  it  can 
be  concluded  that  commercial  polysilicon  cells  are  not  suitable  to  handle  high 
concentration  ratios  unless  cooled  extremely  well.  This  fact  also  sheds  light  on  the 
proprietary  cooling  systems  used  in  many  state-of-the-art  CSP  arrays.  These 
manufacturers  are  using  a  more  effective  cooling  system  than  an  air  fan  to  manage 
thermal  runaway  in  these  cells.  Furthermore,  high-efficiency  solar  cells  are  inherently 
better  at  managing  heat  since  a  higher  percentage  of  the  incident  photon  energy  is 
converted  to  electricity. 
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VIII:  Conclusions 


8.1  Discussion 

This  research  has  demonstrated  the  ability  to  predict  the  steady  state  temperature 
of  a  concentrated  solar  cell.  Furthermore,  the  optimum  concentration  factor  for  the  test 
cell  was  determined  from  simulation.  The  test  results  combined  with  simulated 
predictions  revealed  that  the  optimum  concentration  factor  for  the  test  set-up  was  much 
lower  than  hypothesized,  due  to  the  significant  amount  of  heat  contributed  by  multiple 
suns. 

A  novel  model  was  developed  to  facilitate  the  computer  simulation  of  the  cell 
under  concentration,  and  the  model  was  shown  to  accurately  predict  the  steady  state 
temperature  of  the  cell  at  a  concentration  of  up  to  four  suns.  When  used  to  predict  the 
steady  state  operating  temperature  of  the  cell,  the  Gamma  Model  outperforms  the  other 
electrical  solar  cell  models  discussed  in  Chapter  II  as  shown  in  Figure  37. 


Experimental  and  Simulated  Results  Comparison 


Figure  37:  Electrical  Model  Results  Comparison 
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Although  the  electrical  portion  of  the  model  provides  a  relatively  minor 
component  of  heat  generated  in  the  cell,  using  the  wrong  model  (i.e.  exponential  model) 
can  adversely  affect  the  steady  state  temperature  predictions.  Furthermore,  the 
expression  for  total  current  calculated  using  the  Gamma  Model  captures  the  important 
phenomenon  of  surface  carrier  recombination,  discussed  in  Section  2.8.  This  can  be  seen 
when  the  calculated  photocurrent  is  compared  to  the  current  computed  using  the  ideal 
diode  model  against  the  measured  short  circuit  current  of  the  cell,  shown  in  Figure  38. 


Comparison  of  Photocurrent  Predictions- 1  Sun  Illumination 


Temperature  (K) 

Figure  38:  Simulated  Photocurrent  Comparison 

The  standard  diode  model  photocurrent  Iph  —  qG0A(ln  +  W  +  Lp),  fails  to 
adequately  capture  the  effect  of  surface  recombination  because  the  assumption  is  made 
that  the  majority  of  EHP  generation  takes  place  within  or  near  the  space  charge  region 
[17],  when  in  fact  a  significant  amount  of  carrier  generation  takes  place  in  the  absorption 
layer  of  the  device  [8]. 
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Another  benefit  of  the  Gamma  Model  is  that  the  temperature  dependence  of  the 
photocurrent  is  incorporated  into  the  overall  current  calculation.  Figure  39  demonstrates 
this  contributing  feature  of  using  the  Gamma  Model  in  the  solution  for  the  solar  cell 
current  equation,  showing  that  as  the  temperature  of  the  cell  increases  the  photocurrent 
decreases  because  of  an  increase  in  the  kinetic  energy  within  the  semiconductor  crystal 
lattice.  An  inclusion  of  the  thermal  effects  into  the  electrical  model  yields  a  more 
complete  solution  for  the  dynamic  thermal  effects  of  the  solar  cell  over  the  ideal  diode 
model,  which  assumes  that  the  photocurrent  is  independent  of  changes  in  temperature 
[17]  [12]. 

Temperature  Dependence  of  Gamma  Model  Photocurrent- 1  Sun  Illumination 
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Figure  39:  Temperature  Dependence  Feature  of  Gamma  Model 
Using  the  Gamma  Model  simulation,  an  optimum  concentration  factor  analysis 

reveals  that  an  illumination  level  of  2.9  suns  is  the  best  configuration  for  the  solar  cell 

and  cooling  system  combination  used  in  this  research.  A  higher  concentration  factor 

(more  suns)  contributes  more  heat,  and  as  the  concentration  factor  exceeds  three  suns,  the 
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open  circuit  voltage  begins  to  collapse  and  the  output  efficiency  decreases.  Given  the  test 
cell  and  cooling  system,  it  is  not  practical  to  concentrate  beyond  a  few  suns.  This 
research  has  really  demonstrated  the  importance  of  thermal  management  in  a  CSP 
system. 

8.2  Key  Contributions 

There  are  several  contributions  of  this  research  to  current  and  future  solar  energy 
research  projects.  The  thermal  model  presented  in  this  work  represents  the  first  and  only 
solar  cell  model  that  incorporates  both  the  electrical  characteristics  and  thermal  properties 
of  the  cell  in  order  to  predict  the  steady  state  temperature  under  illumination.  Thennal 
models  developed  from  a  purely  thermodynamic  perspective  were  studied  and  used  in 
this  research,  but  nowhere  in  the  literature  had  someone  combined  the  thermodynamic 
properties  of  the  solar  cell  surface  with  the  solid-state  physics  governing  the  flow  of 
charge  across  the  PN-junction.  This  research,  and  the  Gamma  Model  in  particular,  is 
beneficial  to  future  research  since  different  materials  can  easily  be  incorporated  into  the 
computer  model. 

Another  important  contribution  of  this  research  is  the  novel  model  for  excess 
carrier  concentration  in  a  PN-junction  solar  cell.  Existing  models  use  an  exponential  or 
linear  approximation  on  both  sides  of  the  PN-junction,  but  upon  further  investigation  it 
was  detennined  that  a  more  accurate  model  capturing  the  effect  of  surface  carrier 
recombination  was  necessary  to  complete  this  research.  Furthermore,  it  was  detennined 
that  a  negative  exponential  approximation  of  carrier  concentration  produced  a 
mathematical  singularity  as  the  photon  penetration  depth  approached  the  minority  carrier 


96 


diffusion  length.  Using  a  gamma  function  to  solve  the  second  order  different  transport 
equation  on  the  n-side  yielded  an  initially  complex  solution  that,  when  simplified, 
resulted  in  an  eloquent  solution  to  carrier  concentration.  The  model  was  verified  with 
experimental  results  that  demonstrate  the  ability  of  the  thermal  model  to  predict  steady 
state  temperature  as  a  function  of  concentration.  The  model  can  easily  be  adapted  to 
reflect  different  cooling  systems  and  with  little  manipulation  other  semiconductor 
materials.  However,  the  full  extent  of  these  possibilities  exceeded  the  boundary  of  this 
initial  research. 

8.3  Future  Research 

One  of  the  main  goals  of  this  research  was  to  open  the  door  for  future 
photovoltaic  research  at  AFIT.  As  the  first  modern  solar  power  related  research  project 
at  AFIT,  this  work  has  laid  the  foundation  for  several  future  research  opportunities  in  the 
field  of  solar  energy. 

A  very  exciting  opportunity  lies  in  the  development  of  low-power  methods  to  cool 
CSP  systems.  As  was  shown  in  this  research,  even  a  modest  increase  in  concentration 
causes  a  severe  increase  of  temperature  that  ultimately  degrades  the  output  efficiency  of 
the  cell.  Passive  cooling  methods  such  as  heat  sinks  are  readily  available  to  dissipate 
some  of  the  heat,  but  at  high  concentration  ratios  the  relative  size  of  the  heat  sinks  would 
physically  dwarf  the  cells,  essentially  making  them  infeasible  in  this  capacity  as  power- 
to-weight  ratios  are  considered.  A  real  breakthrough  lies  in  the  discovery  of  an  efficient 
passive  cooling  system  that  requires  no  external  power  and  effectively  manages  the  cell 
temperature  to  an  optimum  operating  condition.  The  goal  of  any  solar  power  system 
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should  be  a  stand-alone  power  supply  that  is  optimized  under  every  illumination 
condition. 

The  test  set-up  used  in  this  research  demonstrated  that  the  active  cooling  system 
used  was  effective  at  cooling  the  cell,  but  only  under  low  concentration  ratios. 
Additionally,  air  was  used  as  the  laminar  fluid  in  this  research,  but  it  would  be  very 
interesting  to  investigate  the  possibility  of  using  a  shallow  liquid  flowing  over  the  surface 
of  the  cell  to  cool  the  system.  A  number  of  transparent  liquid  coolants  exist  to 
accomplish  this  research,  but  again  this  active  cooling  system  would  require  external 
power.  However,  investigation  of  alternative  solar  cell  cooling  methods  (both  active  and 
passive)  would  be  extremely  interesting  and  beneficial  to  the  field  of  solar  energy 
research. 

Another  future  research  opportunity  lies  in  running  the  same  illumination  tests 
with  different  semiconductor  material  solar  cells,  such  as  GaAs,  InGaP,  CdTe,  and  other 
high-efficiency  cell  materials  to  investigate  the  thennal  properties  of  higher  efficiency 
solar  cells,  like  the  ones  used  in  industrial  CSP  systems.  Exotic  multi-junction 
semiconductor  cells  are  used  in  most  modern  CSP  systems,  but  much  of  the  research 
behind  the  design  and  implementation  of  these  systems  is  proprietary. 

A  significant  hurdle  facing  solar  cell  researchers  in  academia  is  that  many 
breakthroughs  in  the  field  are  often  closely  guarded  “trade  secrets”.  Obviously,  if 
someone  develops  a  technology  that  could  dramatically  increase  efficiency  and  decrease 
cost  that  would  create  a  major  competitive  advantage  in  this  extremely  competitive  field. 
Consequently,  the  lucrative  field  of  solar  energy  has  created  a  vacuum  of  academic 
literature,  especially  research  concerning  CSP  systems.  It  is  the  author’s  sincere  belief 
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that  solar  energy  research  is  especially  important  within  the  department  of  defense  as  the 
culture  continues  to  shift  towards  energy  awareness.  As  the  military  searches  for  ways  to 
save  money  and  increase  efficiency  across  all  sectors,  defense-minded  solar  energy 
research  will  only  catalyze  this  effort  to  reduce  demand  and  increase  energy  supply. 
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Appendix  A:  MATLAB  Code 


%  this  m-file  computes  the  total  current  of  the  cell  using  ambipolar 
%  transport  equations  and  the  ideal  diode  law.  The  currents  are 
derived 

%  from  expressions  for  the  excess  minority  carrier  concentrations  on 
either 

%  side  of  the  pn  junction.  The  current  will  change  depending  on  the 
%  illumination,  temperature,  what  else? 
h=6 . 62  6e-34 ; 
c=3e8 ; 

%  IV  curve ! ! ! ! 

PD  =  [1.32  20.96  113.48  182.23  234.43  286.01  355.88  386.80  381.78 
492.18  751.72  822.45  842.26  890.55  1077.07  1162.43  1180.61  1212.72 
1180.43  1253.83  1242.28  1211.01  1244  1299  1273  1276  1277  1292  1284  1262 
1261  1255  1240  1243  1233  1188  1228  1210  1200  1181  973  1173  1152  1133 
974  1110  1086  1070  733  1036  1018  1003  988  860  932  923  914  407  857  843 
835  817  807  793  778  217  163  249  231  255  279  529  496  585  486  448  486  500 
100  116  108  155  139  374  383  424  382  383  323  344  345  284  175  2  30  67  59 
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180  161 
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39 
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57 

53 

50 

31 

28 

24 

15 

2]  ; 

lam 

= 

0  . 

3e-6 : 0 . 

01775e-6 : 

2 . 5e 

-6; 

R  vec 

=  [• 

25  . 

.245 

.24 

.235 

.23 

.225  . 

22  . 

21 

.2 

.19  . 

18  . 

.17 

.16 

.15  . 

.14 

.  13 

. 

12 

.  11 

.1 

.1  . 

1  .1 

.  1  . 

1.1.1 

.1 

.1 

.1 

.1  .1 

.1 

.1 

.1 

.1  .1 

.1 

.1 

.  1 

.1 

.1 

.1 

.  1 

.12 

.14  . 

16  . 

18.2. 

22  . 

24 

.26 

.28 

.3]  ; 

R_vec (54:124)=. 3; 
vv=  3e8./lam; 
PDconv  =  PD.*le-4; 


F=PDconv./ (h.*vv)  ;  %  compute  flux  from  power  density  (124  elements, 

ie  wavelengths) 


PDtg  =  [0.0427  .068  1.33  1.845  2.232  3.521  3.9  4.4  5.3  5.8  6.8  7.1  7.3 
7.9  8.3  8.8  9.7  10.2  11.13  11.5  12  12.5  13.06  13.6  14.5  14.9  15.4  15.9 

16.5  17.7  18.1  18.6  19.12  19.5  19.89  20.02  20.41  21.31  21.7  21.82  22.21 

22.6  22.99  23.63  23.63  23.9  24  24.15  24.15  24  23.8  23.4  23.1  22.34 
21.95  21.7  21.57  21.31  20.92  20.79  20.54  20.15  19.76  19.38  18.99  18.73 

18.6  18.22  17.44  17.06  16.67  15.9  15.5  15.12  14.7  14.7  13.96  13.96 
13.58  13.19  13.19  12.8  12.54  12.16  11.77  11.77  11.51  11.38  11.13  10.74 
10.48  10.35  10.09  9.8  9.45  9.2  8.9  8.8  8.5  8.3  8.03  7.9  7.6  7.4  7.13 
6.8  6.36  6.1  5.9  5.7  5.6  5.3  5.2  4.9  4.8  4.6  4.4  4.2  4  3.8  3.7  3.58  3.4 
3.4]; 

PDtgscale=PDtg . *trapz ( PD) /trapz ( PDtg) ; 


PDtgconv=PDtgscale . *le-4  ; 

Ftg= ( PDtgconv . / (h . *vv) ) ; 

for  T  =300:1:700;  %  temperature  in  Kelvin  (open  circuit  voltage 
decreases  while  Isc  increases  as  T  rises) 

A  =  15  ;  %  cell  area  =  1  cmA2 
q  =  1.619e-19;  %  electron  charge  (C) 
k  =  8.62e-5;  %  Boltzmann's  constant 
kl  =  1.38e-23;  %  Boltzmann's  constant  J/k 
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Vt  =  k*T;  %  thermal  voltage  =  0.0258  V 

ni  =  9 . 38el9* ( (T/300) A2) *exp (-6884/T)  ;  %  source  PV  education  "Intrin. 

Carr.  Cone."  (for  Si) 

ND  =  1.7el6;  %  doping  levels 

NA  =  1 . 7  e 1 6 ; 

Dn  =  1350*Vt/ ( (1+ (81*NA) / (NA+3 . 2el8) ) A ( . 5) ) ;  %  diffusion  coeffecient 
for  electrons  (cmA2/s) 

u_n  =  Dn/ (k*T) ; 

Dp  =  480*Vt/ ( (1+ (350*ND) / (ND+1 . 05el9) ) A ( . 5) ) ;  %  diffusion  coeffecient 
for  holes  (both  temp  dependent) 

u _ p  =  Dp/ (k*T) ; 

pno  =  niA2/ND;  %  initial  minority  carrier  cone  on  the  n-side 

(function  of  temperature) 

npo  =  niA2/NA;  %  initial  minority  carrier  cone  on  the  p-side 

es  =  11 . 7*8 . 85e-12  ;  %  relative  permitivity  for  Si 
vbi  =  (k*T)*  log (NA*ND/ni A2 )  ;  %  built  in  voltage 

xb  =  sqrt ( (2*es* (vbi) /q) * (NA/ND) * (1/ (NA+ND) ) ) ;  %  position  of  the  depl 
region  boundary  on  n-side  (w/  the  center  of  the  depl  region  as  ref) 

xp  =  sqrt ( (2*es* (vbi) /q) * (ND/NA) * (1/ (NA+ND) )) ;  %  position  of  the  depl 
region  boundary  on  the  p-side  (w/  depl  region  center  as  ref) 

tau_p  =  (12/ (1+ (ND/5el6) ) ) *le-6  ;  %  minority  carrier  lifetimes  (s) 

tau_n  =  (12/ (1+ (NA/5el6) ) ) *le-6  ;  %  (s) 

Lp  =  sqrt (Dp*tau_p)  ;  %  diffusion  length  (cm) 

Ln  =  sqrt(Dn*tau  n)  ;  %  (cm)  make  sure  units  work  out 

C  =  3;  %  CONCENTRATION  FACTOR! ! 

Eg  =  1.17  -  4.73e-4*  (TA2) / (T+636)  ; 

Eg  max  =  1.1245; 
vmin=Eg  max*q/h; 

lambdamax=1000000*c/vmin;  %  1.1034  um 

dlambda=0 . 01775;  %  to  get  124  elements 

lambda_count=0 ; 

Itot=zeros (1001, 124) ; 
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Itotl=Itot; 

Itot2=Itot; 

Itot3=Itot; 

Itot4=Itot; 

Itot5=Itot; 

for  lambda=0 .3:0.01775:2.5;  % 

lambda  count=lambda  count+1; 


R=  R  vec (lambda_count) ;  %  reflectivity  for  Si  about  20% 
v= ( c / 1 ambda ) *1000000; 

Ftemp=dlambda*Ftg (lambda_count) . *C; 

11  =  0.5; 

12  =  0.8; 

13  =  1.1; 

if (lambda<ll) 

alpha  =  10A (-6 . 7*lambda+7 . 4) ;  %  how  far  photons  penetrate  into  material 
(1/cm) 
end 

if ( (lambda>ll) && (lambda<12) ) 
alpha  =  10 A (-3 . 3*lambda+5 . 6) ; 
end 

if  ( (lambda>12) && (lambda<13) ) 
alpha  =  10A (-6 . 7*lambda+8 . 4) ; 
end 

if  (lambda>13) 
alpha  =  0  ; 
end 

alpha=alpha;  %  to  go  from  1/m  to  1/cm 
alpha  vec (lambda_count) =alpha; 

%  alpha  will  vary  from  10A6  -  10 

pen_depth=  1. /alpha  vec;  %  penetration  depth  in  cm 
%  hv  is  the  number  of  joules  per  photon 

Go_prime  =  alpha*Ftemp  ;  %  gen  rate  at  surface  photons/cmA3/s 
nthick  =  3e-4  ;  %  thickness  of  n-side  (cm) 

Go  =  Go_prime*exp ( - (nthick-xb) . *alpha) * (1-R)  ;  %  gen  rate  at  edge  of 

depl  boundary 

Go_vec (lambda_count) =Go; 

W=xb+abs (xp)  ;  %  width  of  the  depletion  region 

count  =  0  ;  %  store  each  current  in  a  diff  count  variable 
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if ( lambda<=lambdamax) 

for  Va  =  0:0.001:1;  %  this  is  in  VOLTS 

count  =  count  +  1; 

%if ( (1/ (LpA2 ) ) >alphaA2 ) 

Itotl (count, lambda_count) =A* (q*Dp* ( (pno* (exp (Va/ (k*T) ) - 
1) *exp (2*xb/Lp) ) /  (1  -exp (2*xb/Lp) ) * (1/Lp) ) ) ; 

Itot2 (count, lambda_count) =A* (q*Dp* ( (pno* (exp (Va/ (k*T) ) -1) ) / (1  - 
exp (2*xb/Lp) ) * (1/Lp) ) ) ; 

Itot3 (count, lambda_count) =A*q*Dp* ( (Go . * (1/ alpha+ (nthick- 
xb) )  )  ./  (2*Dp) ) ;  %new  photo  current 

% I tot 3 (count, lambda_count ) =A*q*Dp* (alpha . *Go . / (Dp* ( 1 / (LpA2 ) - 
alphaA2)));  %photo  current  flows  in  the  reverse  bias  direction  (opposi 
Itot4 (count, lambda_count) =-A*q*Dn* (npo* (exp (Va/ (k*T) ) - 
1) *exp (xp/Ln) *1/Ln) ; 

%end 

%Itot5 (count, lambda  count) =A*Go*q* . 01  ;  %  any  photon  that  makes  it  to 
the  depletion  region  gets  turned  into  a  carrier 

Itot5 (count, lambda  count) = (Go_prime . / alpha) *A*q* (exp (-nthick+xb) -exp (- 
nthick-xp) ) ; 

Itot (count, lambda_count) 

=Itotl (count, lambda  count) +Itot2 (count, lambda_count) +Itot3 (count, lambda 
_count) +Itot4 (count, lambda_count) +Itot5 (count, lambda_count) ; 

Vaplot (count)  =  Va  ; 

end 

end 

end 

Itot22=sum ( Itot, 2 )  ; 
hold  on 

P  =  Vaplot .* Itot22  '  ;  % 

xx=f ind ( P>=0 ) ; 
if  (T  ==  361) ; 

plot (Vaplot (xx) , Itot22 (xx) ) ; 
end 

maxpowercount  =  find(P  ==  max (P) ) ; 

Iop(T-299)  =  Itot22 (maxpowercount)  ;  %  operating  current  for  max  power 


%%  This  subroutine  m-file  computes  the  entire  time  dependent  component 
of 

%  heat  generated  by  the  cell  under  illumination.  In  order  to  reach 
steady 

%  state,  the  amount  of  heat  generated  by  the  cell  must  eventually  be 
%  balanced  out  by  the  amount  of  heat  dissipated  by  the  cell 

nO  =  (ND-NA) / 2  +  sqrt ( ( (ND-NA) /2 ) A2+ni A2 )  ;  %  niA2  =  n*p  ;  NA+n  = 

ND+p 

pO  =  niA2/nO  ; 


104 


sigma  =  q* (nO*u  n+pO*u_p)  ;%  mobility  of  both  holes  and  electrons 
DECREASES  with  increasing  temp! 

rho=  1/sigma  ;  %  sigma  is  conductivity,  increases  with  temp  increase 

W=  1  ;  %  cm 
L=  1  ;  %  cm 

pthick  =  .06  ;  %  thickness  in  cm 

t=  470e-6  ;  %  total  cell  thickness 

%Rcell  =  rho*L/ (W*t)  ;  %  ohm-cm  sheet  resistance 

Rcell  =  .10;  %  ohms 

Pr=  lop (T-299) A2*Rcell  ;  %  heat  due  to  illumination  (Itot  computed  in 
another  subroutine) 


Eg  =  1.17  -  4.73e-4* (TA2) / (T+636)  ;  %  source 

http : // www .ioffe.rssi.ru/ SVA/NSM/ Semicond/ Si /bands tr . html 

lambda_count2=0 ; 

%  this  is  for  photons  ABOVE  bg  energy 

for  lambda3  =  0 . 3 : 0 . 01775 : lambdamax; 
lambda_count2=lambda_count2+l ; 

R=R_vec ( lambda  count2); 

Ftemp=dlambda*Ftg ( lambda  count2) .  *C; 


vv 2= (c/lambda3) *1000000; 

Pa2 ( lambda_count2 ) = (A) . * ( ( (h* (vv2) /q) -Eg) *q* (1-R) . * (Ftemp) ) ; 

end 

for  lambda2=lambdamax+ . 01775 : . 01775 : 2 . 5;  %  1.1034  is  lambdamax  (140 
data  pts) 

lambda  count2=lambda  count2+l; 

R=R  vec ( lambda_count2 ) ; 

Ftemp=dlambda*Ftg ( lambda_count2 ) . *C; 

vv= (c/lambda2) *1000000; 

%Pa2 (lambda  count2)=  (A) * (h*c/lambda2) * (Ftemp) * (1-R)  ; 

%  accounts  for  temp  change  F=  flux  (J/sec) !  increases  as  temp  rises 
since  the  band  gap  decreases 

Pa2 ( lambda_count2 ) =  (A) . * ( ( (h* (vv) /q) ) *q* (1-R) . * (Ftemp) ) ; 
end 

%Pa_tot=sum ( Pa2 ) +sum (Pa) ; 

Pa_tot=sum ( Pa2 )  ; 

Pgen  =  Pr  +  Pa  tot  ; 

%Pgen=  Pa_tot; 

%%  19  July  2011 

%  This  subroutine  m-file  computes  the  time  dependent  amount  of  heat 
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%  dissipated  from  the  solar  cell  under  illumination 

sigma  =  5.6704e-8  ;  %  W / (mA2*KA4) 

Pbb  =  (A) * . 0001*sigma* (TA4) ;  %  blackbody  radiation  (convert  area  to 
mA2) 

Tstream  =  300  ; 

Hsil=0 .0123; 

%tau_H=9 . 259;  %  1.75  m/s  (slowest  setting  on  the  POT) 
tau_H=6.625;  %  4.25  m/s  (fastest  setting  on  POT) 

%tau  H=l;  %  example 
hplate  =  Hsil/tau  H  ; 

Pc  =  hplate* (A) * (T-Tstream)  ;%  convective  cooling  of  the  cell  by  the 
ambient  air 

%  as  it  stands  now.  Pc  remains  constant  with  temp  (could  be  refined  to 
%  reflect  changes  in  T) 

Pdiss  =  Pbb  +  Pc  ;  %  total  heat  dissipated  from  the  solar  cell 

o,  o, 
o  o 

Ptot(T-299)  =  Pgen  -  Pdiss  ; 
end 

Ptot2=abs (Ptot) ; 

steadystatecount  =  find(Ptot2  ==  min(Ptot2)); 
steady_state_temp=steadystatecount+2  99 
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Appendix  B:  Convective  Cooling  Coefficient  Derivation 


The  expression  for  Hpiate  shown  by  equation  (6-1)  was  derived  using 
thennodynamic  steady  state  cooling  properties.  Initially,  the  change  in  temperature  with 

respect  to  time  ^  is  expressed  in  terms  of  the  total  power  dissipated  and  the  inverse  of 

specific  heat,  where  Hsif  is  the  inverse  of  specific  heat  (for  Si)  and  PD  is  the  total  power 
dissipated. 


dT 

~dt  =  Hsil  PD 


(B-1) 


The  resulting  dimensional  analysis  of  (B-1)  proves  the  validity  of  this  expression, 
where  temperature  is  in  units  of  Kelvin  K,  time  is  in  units  of  seconds,  power  is  in  units  of 
Joules  per  second,  and  the  inverse  of  specific  heat  is  in  units  of  Kelvin  per  Joule  . 


(B-2) 


As  demonstrated  in  Chapter  IV,  specifically  equations  (4-11)  and  (4-12),  the  heat 
dissipated  can  be  expressed  as  a  function  of  both  the  convective  cooling  and  the 
blackbody  cooling.  Additionally,  the  change  in  temperature  with  respect  to  time  should 
be  negative  since  the  convective  cooling  will  decrease  the  temperature  with  time.  Using 
the  full  expression  for  power  dissipated  and  the  correct  sign  for  the  change  in 
temperature,  equation  (B-1)  can  be  more  completely  expressed  by  equation  (B-3). 


~dt=  n^^viateV  -  Tc)  +  o{T  -  Tcy) 

Equation  (B-3)  takes  the  form  of  an  ordinary  differential  equation,  and  must  be 
solved  in  order  to  extract  an  expression  for  temperature.  Following  the  procedure  used  to 
solve  any  differential  equation,  a  conjecture  to  the  solution  of  equation  (B-3)  is  made 
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using  an  exponential  form  for  temperature,  given  by  equation  (B-4)  where  Tm  is  the 
maximum  temperature,  Tc  is  the  temperature  of  the  convective  cooling  stream,  and  rH  is 
the  exponential  cooling  decay  rate. 


_ t_ 

T  =  (Tm-  Tc)e  +  Tc 


(B-4) 


Referring  to  equation  (B-4),  at  t  =  0  the  light  source  is  turned  off,  and  at  t  =  oo 
the  steady  state  temperature  point  is  reached  and  the  temperature  is  equal  to  Tc.  The 
solution  conjecture  given  by  (B-4)  must  be  used  in  the  original  differential  equation  to 
detennine  its  validity  as  a  solution. 


e 


_ £_  _4 1_ 

=  Hsu  Hpiate(Tm  —  Tc)e  Th  +  Hsn  &(Tm  ~~  fc)46  Th 


(B-5) 


Expression  (B-5)  can  be  simplified  by  removing  the  like  exponential  terms  and 
the  like  temperature  difference  terms  (Tm  -  Tc). 


~  Hsu  Hpiate  +  Hsn  <j(Tm  Tc)3e  th 

Since  the  goal  of  this  analysis  is  to  detennine  steady  state  temperature,  equation 
(B-6)  can  be  used  to  solve  for  Hpiate  as  t  oo,  since  the  exponential  term  is  driven  to 
zero.  The  resulting  solution  concludes  the  derivation. 


(B-7) 

(B-8) 
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Appendix  C:  Thermal  Coefficient  for  Different  Air  Velocities 


Convective  Cooling  Calibration  Test  #1 


Convective  Cooling  Calibration  Test  #2 
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Convective  Cooling  Calibration  Test  #3 


Time  (sec) 


Figure  42:  Cooling  Coefficient  Test-  Air  Velocity  3.0  m/s 


Convective  Cooling  Calibration  Test  #4 
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Convective  Cooling  Calibration  Test  #5 


Figure  44:  Cooling  Coefficient  Test-  Air  Velocity  1.75  m/s 


-3  H  .  .  as  a  Function  of  Air  Speed 

x  10  P'ate 
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